For  Reference 


NOT  TO  BE  TAKEN  FROM  THIS  ROOM 


6Cx 


m~m 

mmmmm 


The  University  of  Alberta 
Printing  Department 
Edmonton,  Alberta 


Digitized  by  the  Internet  Archive 
in  2020  with  funding  from 
University  of  Alberta  Libraries 


https://archive.org/details/RRao1969 


/ 


THE  UNIVERSITY  OF  ALBERTA 


STABILITY  STUDIES  ON  NONLINEAR  DISCRETE-DATA  SYSTEMS 
THROUGH  FUNCTIONAL  ANALYSIS 


by 


RAMAPPA  SRIDHARA  RAO 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  DOCTOR  OF  PHILOSOPHY 


DEPARTMENT  OF  ELECTRICAL  ENGINEERING 


EDMONTON,  ALBERTA 


FALL,  1969 


3) 


THE  UNIVERSITY  OF  ALBERTA 

FACULTY  OF  GRADUATE  STUDIES 


The  undersigned  certify  that  they  have  read, 
and  recommend  to  the  Faculty  of  Graduate  studies  for 
acceptance,  a  thesis  entitled  STABILITY  STUDIES  ON 
NONLINEAR  DISCRETE-DATA  SYSTEMS  THROUGH  FUNCTIONAL 
ANALYSIS  submitted  by  RAMAPPA  SRIDHARA  RAO  in  partial 
fulfilment  of  the  requirements  for  the  degree  of 


Doctor  of  Philosophy. 


(i) 


ABSTRACT 


This  thesis  deals  with  the  study  of  bounded-input  bounded- 
output  (BIBO)  stability  of  certain  classes  of  nonlinear  discrete- 
data  systems  using  the  methods  of  Functional  Analysis.  Sufficient 
conditions  for  the  BIBO  stability  of  such  systems  are  established 
employing  the  contraction  mapping  principle.  The  systems  considered 
include  those  with  the  nonlinearity  in  the  feedback  path,  as  well  as 
in  the  forward  path.  A  system  with  a  slope  restricted  nonlinearity 
is  also  studied.  The  method  for  finding  the  input  and  output  bounds 
for  stability  in  a  special  case  of  nonlinearity  is  discussed.  Finally 
a  method  for  generating  the  system  solution  by  a  process  of  successive 
approximations  is  given.  It  is  shown  that  the  solution  comes  out  in 
the  form  of  a  discrete  Volterra  type  series,  whose  convergence  is 
assured  in  the  region  in  which  the  BIBO  stability  conditions  derived, 


are  satisfied. 


(ii) 


ACKNOWLEDGEMENTS 


The  author  is  deeply  indebted  to  Dr.  G.S.  Christensen,  the 
supervisor  of  this  project,  for  his  invaluable  assistance  and  guidance 
throughout  the  course  of  this  work.  He  would  like  to  thank 
Dr.  J.W.  Macki,  of  the  Department  of  Mathematics  for  constructive 
criticisms  and  helpful  suggestions,  and  Professsor  Y.J.  Kingma, 
for  initiating  him  into  the  study  of  sampled-data  systems. 

Further  he  wishes  to  acknowledge  with  gratitude  the  help 
and  encouragement  he  received  from  Dr.  G.B.  Walker,  the  Head  of 
this  Department,  throughout  his  stay  in  this  University.  Special 
thanks  go  to  Dr.  V.  Gourishanker ,  who  has  been  largely  responsible 
for  the  author’s  decision  to  launch  on  this  academic  adventure, 
and  whose  friendly  counsel  has  always  been  a  source  of  great 
encouragement . 

Acknowledgement  is  gratefully  given  to  the  National 
Research  Council  of  Canada  for  providing  assistance  for  the  summer 
of  1967,  and  to  the  University  of  Alberta  for  awarding  the  Teaching 
Assis tantship  (1967-68),  the  Intersession  Bursary  (1968  Summer)  and 
the  Dissertation  Fellowship  (1968-69). 


~v 

i 


(ill) 


TABLE  OF  CONTENTS 

Page 

LIST  OF  TABLES .  (v) 

LIST  OF  ILLUSTRATIONS .  (vi) 

LIST  OF  MAJOR  SYMBOLS .  (vii) 

CHAPTER  I:  INTRODUCTION  .  1 

CHAPTER  II:  STABILITY  OF  LINEAR  DISCRETE-DATA  SYSTEMS  ...  5 

2.1  Introduction  .  5 

2.2  Conventional  Time  Domain  and  Z-Domain 
Representation  of  Linear  Discrete-Data 

Systems  .  5 

2.3  The  Functional  Representation  .  10 

2.4  Stability  of  Linear  Discrete-Data 

Systems  .  13 

2.5  Computation  of  ||  B*  ||  for  a  given  Linear 

Discrete-Data  System  .  16 

CHAPTER  III:  STABILITY  OF  NONLINEAR  DISCRETE-DATA  SYSTEMS.  .  21 

3.1  Introduction .  21 

3.2  System  A:  Description  and  Assumptions  .  .  21 

3.3  Functional  Representation  of  System  A  .  .  23 

3.4  Conditions  for  the  BIBO  Stability 

of  System  A .  26 

3.5  Conditions  for  the  BIBO  Stability  of 

System  B  .  33 

CHAPTER  IV:  COMPUTATION  OF  INPUT  AND  OUTPUT  BOUNDS  ....  39 

4.1  Introduction .  39 

4.2  Input  and  Output  Bounds  for  System  A 

with  a  Single  n^  Degree  Term  in  the 
Nonlinearity  .  39 

4.3  Input  and  Output  Bounds  for  System  B 

with  a  Single  n^  Degree  Term  in  the 
Nonlinearity  .  47 


,  IS  9 


(iv) 

Page 

CHAPTER  V:  STABILITY  OF  A  DISCRETE-DATA  SYSTEM  WITH  A 

SLOPE  RESTRICTED  NONLINEARITY  .  51 

5.1  Introduction . 51 

5.2  System  Description  .  51 

5.3  Stability  Criterion  for  N  .  52 

5.4  Numerical  Example  1 .  58 

5.5  Comparison  with  some  Existing  Results  .  .  65 

CHAPTER  VI:  GENERATION  OF  SYSTEM  SOLUTION  .  69 

6.1  Introduction .  69 

6.2  Derivation  of  the  Solution .  69 

CHAPTER  VII:  CONCLUSION  .  78 

7.1  Summary .  78 

7.2  Comments  and  Suggestions  .  79 

REFERENCES .  81 


f 


(v) 


LIST  OF  TABLES 

Table  Page 

1  Best  Bounds  for  ^ for  BIBO  Stability 

de 

using  Contraction  Mapping  .  65 

2.  Best  Bounds  for  for  BIBO  Stability 

de 

(Jury’s  Method  vs  Contraction  Mapping)  .  .  67 


£  '  >  A 


(Vi) 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

2.1  A  Linear  Discrete-Data  System  .....  .  5 

3.1  Nonlinear  Discrete-Data  System  A  .  . .  22 

3.2  Linearized  System  Corresponding  to  System  A  ...  23 

3.3  Nonlinear  Discrete-Data  System  B  .  33 

3.4  Linearized  System  Corresponding  to  System  B  ...  33 

4.1  Discrete-Data  System  with  a  Cubic  Nonlinearity 

in  the  Feedback  Path .  44 

4.2  Discrete-Data  System  with  a  Cubic  Nonlinearity 

in  the  Forward  Path .  49 

5.1  Class  N  Nonlinear  Discrete-Data  System 

(a)  System  Configuration  (b)  Gain  Function  .  .  51 

5.2  Nonlinear  Integral-Control  Sampled-Data  System  .  .  58 

6.1  Discrete-Data  System  with  a  Cubic  Nonlinearity  in 

the  Feedback  Path .  69 


-\ 

(vii) 


LIST  OF  MAJOR  SYMBOLS 


r 

II 

rt 

'w' 

** 

continuous  input. 

c 

II 

o 

rt 

continuous  output. 

e 

=  e(t)  , 

continuous  error. 

r* 

=  r*(t)  , 

sampled  input. 

c* 

=  c*(t)  , 

sampled  output. 

e* 

=  e*(t)  , 

sampled  error. 

T  Period  of  sampling. 

x(nT)  Value  of  sampled  variable  x  at  the 

sampling  instant  nT  ,  n  =  0,  1,  . .  . ,  . 

6 (t-nT)  Unit  impulse  occurring  at  t  =  nT  . 

=  1  when  t  =  nT,  and  0  when  t  ^  nT. 

00 

x*  =  x*(t)  =  1  x(nT)6(t  -  nT)  ,  using  the  6  notation. 

n=o 


s 

x(s) 

x*(s) 

z 

x(z) 

G(s) 


=  {x  }  ,  where  x  =  x(nT)  ,  using  the  sequence 

n  n=o  n 

notation . 

Laplace  Transform  variable. 

Laplace  Transform  of  the  continuous  function  x(t)  . 
Laplace  Transform  of  the  sampled  function  x*(t) 

Z- transform  variable. 

Z-transform  of  x*(s)  . 

Transfer  function  of  the  linear  plant  G  . 


*A 

J 


(viii) 


G(z) 

Z-transfer  function  of  the  linear  plant  G  . 

H(z) 

Overall  Z-transfer  function  of  the  linearized  system. 

h(kT) 

Value  of  the  sampled  impulse  response  of  the  linearized 

system  at  t  =  kT  ,  k  =  0 ,  1,  . .  .  ,  00  . 

d(x1,x2) 

Distance  between  x^  and  . 

F 

1  *  1 

1 1  *  1 ! 

Banach  space. 

Absolute  value  sign. 

Norm  sign. 

R* 

=  ||r*||,  norm  of  r*  . 

c* 

=  |j  c*  || ,  norm  of  c*  . 

B* 

Convolution  summation  operator. 

H* 

=  jB*  ,  norm  of  the  operator  B*  . 

A 

Nonlinear  operator. 

X 

Lipchitz  constant,  0  <  A  <  1  . 

a 

radius  of  the  sphere  S  . 

e 

Belongs  to. 

< 

Less  than 

< 

Less  than  or  equal  to. 

> 

Greater  then 

> 

Greater  than  or  equal  to 

Sgn 

_ .  _  X 

Signum,  Sgn  x  =  . 

1  x  1 

BIBO 

Bounded-input  bounded-output. 

e 

Arbitrarily  small  real  number  >  0  . 

CHAPTER  I 


INTRODUCTION 


Discrete-data  systems,  or  sampled-data  systems  as  they  are 
commonly  called,  have  come  to  play  a  significant  role"  in  present-  day 
control  technology.  Although  the  use  of  sampled-data  in  control  systems 
was  recognized  quite  early  in  the  history  of  feed-back  control  systems, 
interest  in  the  analysis  and  design  of  sampled-data  systems  began  perhaps 
in  the  early  forties,  when  during  the  war,  problems  connected  with  radar 
tracking  systems  were  encountered.  The  advent  of  the  digital  computers 
and  their  incorporation  into  control  systems  have  given  a  fillip  to 
the  study  of  sampled-data  systems,  as  the  basic  operation  of  such 
computers  by  necessity  tends  to  reduce  them  to  sampled-data  systems. 

An  important  area  of  investigation  into  the  behavior  of  these 
systems  is  the  determination  of  conditions  for  system  stability.  For 
this  purpose  sampled-data  systems  may  be  broadly  divided  into  two 
categories:  the  linear  sampled-data  systems  and  the  nonlinear  sampled- 

data  systems.  For  linear  sampled-data  systems  the  formulation  of  both 
necessary  and  sufficient  conditions  for  system  stability  have  been 
exhaustively  investigated.  However,  many  sampled-data  systems  of 
practical  importance  are  nonlinear  for  which  formulation  of  conditions 
for  system  stability  is  far  from  being  complete,  although  it  has  focussed 
the  attention  of  many  recent  investigators. 

Much  of  the  work  done  in  this  field  deals  with  the  direct 
application  of  Liapunov's  second  method  extended  to  systems  of  difference 
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equations  .  Although  this  is  the  most  general  of  all  known  methods  of 
stability  analysis  of  nonlinear  systems,  it  suffers  from  some  arbitrariness 
in  the  choice  of  suitable  functions. 

A  few  years  ago  the  Rumanian  scientist  V.M.  Popov  developed 
an  entirely  new  approach  to  the  classical  problem  of  absolute  stability 
by  working  in  the  frequency  domain.  Popov’s  results  which  originally 
applied  to  autonomous  continuous  systems  with  nonlinearities  confined 

to  specified  gain  sectors,  have  been  extended  to  discrete-data  systems 

2  3-5 

by  Tsypkin  .  Jury  and  Lee  further  extended  the  application  of 
Popov’s  method  to  nonlinear  sampled-data  systems,  and  by  placing 
constraints  on  the  slope  of  the  nonlinearity  obtained  less  conservative 
results  than  did  Tsypkin.  All  these  stability  criteria  could  only 
predict  absolute  stability  of  the  null  solution,  which  means  that 
sampled  inputs  which  tend  to  zero  fast  enough  produce  outputs  which 
tend  to  zero. 

Another  very  important  kind  of  stability  is  the  requirement 
that  bounded  inputs  produce  bounded  outputs,  usually  referred  to  as 
bounded- input  bounded-output  (BIBO)  stability.  This  thesis  is  primarily 
concerned  with  the  stability  of  certain  classes  of  nonlinear  discrete- 
data  systems  in  the  BIBO  sense.  For  linear  systems  uniform  asymptotic 
stability  of  the  autonomous  system  implies  BIBO  stability  as  a  direct 
consequence.  In  general  this  is  not  true  for  nonlinear  systems  and 
many  counterexamples  have  been  constructed  .  Also  a  nonlinear  system 
can  exhibit  local  BIBO  stability  even  though  it  may  not  be  BIBO 
stable  in  a  global  sense. 

* 


Numbers  placed  above  the  line  of  text  refer  to  the  references. 
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In  a  recent  paper  Iwens  and  Bergen  have  shown  that  the 

3 

condition  given  by  Jury  and  Lee  for  the  absolute  asymptotic  stability 
of  certain  classes  of  autonomous  nonlinear  sampled-data  systems  establishes 
absolute  BIBO  stability  as  well.  Apart  from  this  reference,  available 
technical  literature  provides  little  evidence  of  much  work  having  been 
done  in  the  area  of  BIBO  stability  of  nonlinear  discrete-data  systems. 


In  this  thesis  the  problem  is  approached  in  an  altogether 

different  way.  The  method  of  analysis  is  based  on  the  fixed-point 

property  of  certain  contraction  mappings  in  Banach  space.  Although  this 

concept  is  well  established  in  Functional  Analysis  and  cognate  branches 

of  modern  mathematics,  its  application  to  the  study  of  stability  and 

convergence  problems  arising  in  nonlinear  system  analysis  is  relatively 

new.  Observations  on  the  use  of  this  approach  has  been  made  by  several 

18  9 

authors  in  the  past  -  Kalman  and  Bertram  ,  Zames  ,  and  Sandberg  ,  to  name 

a  few.  On  the  application  side  papers  by  Desoer^,  Leon  and  Ande r son 
12 

and  Leibovic  may  be  mentioned.  Recently  the  BIBO  stability  of  a  class 

of  nonlinear  continuous  data  systems  has  been  successfully  investigated 
13 

by  Christensen  using  the  same  principle.  However,  as  far  as  the  BIBO 

stability  of  nonlinear  discrete-data  systems  is  concerned,  this  approach 

does  not  seem  to  have  been  attempted  so  far,  though  methods  based  on 

contraction  mapping  have  been  used  in  a  different  way  to  investigate 

14  15 

asymptotic  stability  in  the  large  of  such  systems  * 


The  objective  of  this  thesis  is  to  extend  the  technique  of 
contraction  mapping  to  obtain  sufficient  conditions  for  the  BIBO 
stability  of  different  classes  of  nonlinear  discrete-data  systems. 
Some  of  the  results  have  been  presented  by  the  author  in  a  paper 


■ 
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entitled  "Contraction  mapping  applied  to  stability  analysis  of  a  class 

of  nonlinear  discrete-data  systems",  at  the  Second  Asilomar  Conference 

16 

on  Circuits  and  Systems,  October  1968 

Chapter  II  considers  the  stability  of  linear  discrete-data 
systems  from  a  functional  point  of  view.  The  mathematical  concepts 
and  functional  nontations  involved  are  introduced  in  this  chapter. 

Chapter  III  deals  with  the  stability  analysis  of  nonlinear 
discrete-data  systems  with  a  polynomial  type  of  nonlinearity.  The 
contraction  mapping  theorem  is  stated  and  proved,  and  its  application 
for  obtaining  sufficient  conditions  for  the  BIBO  stability  of  the 
systems  considered  is  explained. 

In  Chapter  IV  the  method  of  computation  of  input  and  output 
bounds  in  a  special  case  of  the  nonlinearity,  is  discussed.  Numerical 
examples  illustrating  the  procedure  are  also  given. 

The  analysis  of  systems  with  a  slope  restricted  nonlinearity 

is  treated  in  Chapter  V.  A  criterion  for  the  BIBO  stability  of  such 

a  system  is  stated  and  proved.  Its  application  is  illustrated  with  the 

aid  of  a  numerical  example  and  the  results  compared  with  those  obtained 

3 

by  the  criterion  of  Jury  and  Lee  . 

Chapter  VI  considers  the  generation  of  the  actual  system 
solution  by  a  process  of  successive  approximations,  leading  to  a  discrete 
Volterra  type  series,  and  shows  how  the  stability  conditions  derived 
establishes  the  convergence  of  such  a  series. 

In  conclusion  the  advantages  and  disadvantages  of  the  method 
of  contraction  mapping  are  discussed,  and  suggestions  for  future 
investigations  in  this  field  given. 


CHAPTER  II 


STABILITY  OF  LINEAR  DISCRETE-DATA  SYSTEMS 


2 . 1  Introduction. 

The  main  aim  of  this  chapter  is  to  discuss  the  BIBO  stability 
of  linear  discrete-data  systems  from  a  functional  point  of  view,  as  this 
forms  the  basis  for  the  development  of  the  conditions  for  the  stability 
of  the  nonlinear  discrete-data  systems  considered  later.  This  is 
essentially  a  time  domain  approach,  and  involves  the  representation  of 
the  system  as  a  linear  mapping  in  a  function  space  and  defining  a  norm 
for  the  linear  operator  which  effects  the  transformation.  Before 
introducing  the  functional  notation,  however,  the  more  conventional 
time  domain  representation  of  the  system  will  be  first  discussed,  as 
also  the  Z-domain  representation. 


2 . 2  Conventional  Time  Domain  and  Z-Domaln  Representation  of  Linear 

17-19 

Discrete-Data  Systems 


A  conventional  linear  single-loop  discrete-data  system  is 
shown  in  Fig.  2.1. 


c*(t) 

c(t) 


Fig.  2.1  A  Linear  Discrete-Data  System. 
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G(s)  is  the  transfer  function  of  the  linear  plant  assumed 
to  be  situated  in  the  forward  path.  The  sampler  with  uniform  period  T  , 
and  shown  as  a  switch,  is  a  device  such  that: 


e*(t)  =  e(t)  J  6(t-nT) 
n=o 


l  e(nT)6(t-nT)  . 
n=o 


(2.1) 


That  is  it  produces  a  train  of  impulses,  the  area  of  which  being  the 
value  of  the  sampled  variable  at  the  sampling  instant  nT  .  Note  that 
the  input  to  the  sampler  is  the  continuous  error  function  e(t)  assumed 
to  be  zero  for  t  <  0  ,  whereas  the  output  is  a  discrete  valued  function 
of  t  ,  defined  at  the  discrete  instants  t  =  nT,  n  =  0,1, 2,..., 00 

If  now  equation  (2.1)  is  Laplace  transformed  one  gets: 


e*(s)  =  j  e(nT)e 


-nTs 


L 

n=o 


(2.2) 


since  the  Laplace  transform  of  the  impulse  6(t-nT)  is  equal  to  e 

Ts 

Letting  e  =  z  ,  equation  (2,2)  may  be  written  as: 


•nTs 


e  *  (  s )  I. 


=  e  (  z )  = 


s  =  y  lo§  z 


CG 

'u  e(nT)z  n 
n=o 


(2.3) 


Equation  (2.3)  defines  the  Z-transform  of  the  input  function  e(t)  . 

Consider  now  the  forward  path  of  the  system  shown  in  Fig.  2.1. 
The  input  to  G(s)  is  e*(t)  ,  a  train  of  impulses  with  amplitudes 
e(nT)  at  t  =  nT  ,  whereas  the  output  is  c(t)  ,  a  continuous 
function  of  time.  Since  our  analysis  is  confined  to  the  discrete 


4 


7 


function  space,  a  fictitious  sampler  is  introduced  at  the  output 
synchronized  with  the  input  sampler  for  obtaining  the  sampled  output, 
which  is  denoted  by  c*(t)  .  A  relation  between  c*(t)  and  e*(t) 
may  be  obtained  as  follows: 

Suppose  a  unit  impulse  is  applied  to  G(s)  at  t  =  0  ,  the 
output  response  would  be  g(t)  ,  the  impulse  response  of  G  ,  and  the 
output  of  the  fictitious  sampler  would  be: 

oc 

g*(t)  =  l  g (nT) 6 ( t-nT)  .  (2.4) 

n=o 

Hence  if  e*(t)  is  applied  to  G  ,  the  output  sample  c(nT)  at 
t  =  nT  ,  would  be  the  effects  of  all  samples  e(nT)  ,  e(nT-T)  , 

...,  e(T),  e(0)  .  That  is: 

c  (nT)  =  e(0)g  (nT)  +  e(T)g(nT-T)  +  ...  +  e(nT)g(0)  .  (2.5) 

Equation  (2.5)  may  also  be  written  as: 


n 

c (nT)  =  l  e (kT)g (nT-kT) 
k=o 


n 

=  l  g(kT)e (nT-kT) 
k=o 


c*(t)  will  be  then  given  by: 


c*(t) 


I 

n=o 


n 

I 

L  k=o 


I 

n=o 


n 

l 

L  k=o 


e(kT)g (nT-kT) 


g(kT)e (nT-kT) 


6 (t-nT) 


6 (t-nT) 


(2.6) 


(2.7a) 


(2.7b) 


Thus  the  output  samples  can  be  obtained  from  the  input  samples 


, 
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by  a  process  of  summation  defined  by  equation  (2  0  7a)  or  ( 2  0  7b ) . 
operation  is  called  the  convolution  summation. 

The  foregoing  can  be  related  to  the  Z-domain  approach 
following  way: 

Taking  the  Laplace  transform  of  equation  (2,4)  yields: 


G*(s) 


l  g (nT) e  nTs 
n=o 


l .  e . 


G(z)  =  l  g (nT) z  n 
n=o 


Equation  (2.9)  defines  the  Z-transfer  function  of  G(s)  . 


Now  multiplying  both  sides  of  equation  (2.5)  by  e 


-nTs 


and  taking  the  summation  for  n  =  0  to 


co  cx>  oo 

l  c(nT)e“nTs  =  l  e(0)g(nT)e"nTs  +  £  e(T)g(nT-T) 


-nTs 


n=o 


n=o 


n=o 


+  . . .  +  l  e (nT) g (0)e 
n=o 


-nTs 


or 


l  c(nT)e 
n=o 


nTs  =  l  g (nT) e  nTs  [ e ( 0 )  +  e(T)e"Ts  + 
n=o 

-2Ts 

+  e(2T)e  +  . . 0] 


from  which, 

00  00  oc 

l  c (nT) e_nTs  =  l  g(nT)e“nTs  l  e(nT)e"nTs 

n=o  n=o  n-o 


or 


This 

in  the 

(2.8) 

(2.9) 


(2.10) 

(2  o 11) 

(2.12) 
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c*(s)  -  G* (s)e* (s)  . 

Using  the  Z-transform  notation,  (2.13)  may  be  written  as: 


(2.13) 


c(z)  =  G(z)e(z) 


(2.14) 


Equation  (2.14)  is  the  Z-domain  equivalent  of  equation  (2.7a).  Thus  it 
is  seen  that  convolution  in  the  time  domain  is  equivalent  to  multiplication 
in  the  Z-domain. 

Now  consider  the  closed-loop  system  of  Fig.  2.1.  To  obtain 
the  sampled  response  to  any  arbitrary  input  r(t)  ,  one  may  proceed  as 
follows:  It  is  easy  to  see  that, 

e(z)  =  r (z)  -  c (z)  .  (2.15) 


Substituting  (2.15) 


Puttlng  rv gw 


into  (2„14)  and  simplifying: 

c(z)  =  TTirfe-  ‘  r(z)  • 

=  H(z)  one  may  write: 
c(z)  =  H(z)r(z) 


(2.16) 


(2.17) 


a  relation  which  is  analogous  to  equation  (2C14).  Here  H(z)  is  called 
the  overall  Z-transfer  function  of  the  closed-loop  system.  The  correspond¬ 
ing  time  domain  relationship  between  c*(t)  and  r*(t)  will  then  be 
given  by: 


c*(t)  =  l 
n=o 


n 

l  h (kT) r (nT-kT) 
k^o 


6 (t-nT) 


where  h(kT)  =  Z  ^H(z) 


the  inverse  Z-transform  of  H(z) 


(2.18) 

Note  that 
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h(kT)  is  the  same  as  the  sampled  impulse  response  of  the  closed-loop 
system  at  t  =  kT  .  Given  H(z)  one  method  of  obtaining  h(kT)  is  to 
apply  Cauchy's  integral  formula: 

h  (kT)  =  -r^T  i  H(z)zk_1dz 

ZTTJ  *  _ 

J  p 

k-1 

=  Sum  of  the  residues  of  H(z)z 

at  the  poles  of  H(z)zk  ^  .  (2.19) 

20  21 

2.3  The  Functional  Representation  * 

The  natural  representation  for  any  system  is  a  correspondence 
between  input  and  output  elements.  The  classical  time  domain  relation¬ 
ship  between  the  sampled  input  and  the  sampled  output  given  by  equation 
(2,18)  for  the  linear  discrete-data  system  of  Fig.  2.1,  for  example  is 
one  such  representation.  Here  the  emphasis  is  on  the  response  of  a 
particular  system  to  a  particular  input.  In  order  to  obtain  a  qualitative 
understanding  of  the  properties  and  behavior  of  systems  in  general,  a 
characterization  that  yields  a  mathematical  model  compatible  with  a  large 
class  of  systems  and  a  large  class  of  inputs  is  called  for.  This  is  the 
aim  of  the  functional  approach.  This  approach  utilizes  as  its  tools  the 
concept  of  abstract  spaces  and  operations  on  such  spaces.  The  input 
and  output  functions  are  identified  as  belonging  to  an  abstract  space, 
endowed  with  a  mathematical  structure  suited  to  the  type  of  analysis 
desired,  and  the  behavior  of  the  system  studied  in  this  setting. 

The  space  best  suited  for  the  purpose  of  our  analysis  is  the 
Banach  space.  This  space  is  chosen  as  it  facilitates  the  application 
of  the  contraction  mapping  theorem,  when  the  stability  of  the  nonlinear 


? 
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discrete-data  systems  is  considered.  The  following  definitions  will 
clarify  what  is  meant  by  a  Banach  space. 

Definition  1:  A  collection  of  elements  together  with  a  certain 
structure  of  relations  between  elements  or  of  rules  of  manipulation  and 
combination,  the  whole  supporting  a  mathematical  development  is  often 
called  a  space. 

Definition  2:  A  space  X  is  said  to  be  a  linear  (or  vector) 
space  if  addition  and  scalar  multiplication  are  defined  on  X  satisfying 
the  commutative,  associative  and  distributive  laws.  The  scalar  multi¬ 
plication  is  related  to  some  associated  field  of  scalars,  usually  denoted 
by  . 

Definition  3:  A  space  X  is  said  to  be  a  metric  space,  if 
we  associate  with  any  two  points  x^ ,  x ^  of  X  a  distance  function  or 
metric  such  that: 


(a)  d(x  ,x2)  >  0  ,  if  x1  ¥  x2  . 

(b)  d(x^,x2)  =  0  ,  if  and  only  if  x^  =  x^  . 

(c)  d(xisx  )  =  d(x  ,x  )  . 

(d)  d(x1,x2)  <_  d(x^ >x^)  +  d (x ^ , x 2 )  for  every  x3  e  X  . 


Definition  4:  A  norm  on  a  linear  space  X  is  a  real  valued 
function,  whose  value  at  x  we  denote  by  ||x||  ,  with  the  properties: 

(a)  ||  x  ||  _>  0  . 

(b)  || x ||  ¥  0  ,  if  x  i  0  . 

(c)  ||x1  +  x2  ||  <_  1 1 x x  j|  +  1 1 x 2  ||  . 

(d)  |jaxj|  =  |  a  j  [j  x  1 1 ,  for  all  a  e  u?" . 
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Definition  5:  A  linear  space  on  which  a  norm  is  defined  becomes 
a  metric  space  if  we  define  dCx^jX^)  =  |jx  -  x  ||  .  A  linear  space  which 
is  a  metric  space  in  this  way  is  called  a  normed  linear  space  or  a  normed 
vector  space. 

Definition  6:  A  sequence  {x^}  in  a  metric  space  is  called 
a  Cauchy  sequence  if  d(xn,x^)  0  as  m  and  n  -*■  00  . 

Definition  7:  If  a  metric  space  X  has  the  property  that 
every  Cauchy  sequence  in  X  has  a  limit  in  X  ,  X  is  said  to  be 
complete. 

Definition  8:  If  a  normed  linear  space  is  complete  it  is 
called  a  Banach  space. 

The  input  and  output  functions  r*(t)  and  c*(t)  are  real 
valued  functions  of  time  defined  only  at  the  discrete  instants  t  =  nT  , 
n  =  1,2,..., 00  ,  on  the  interval  <  t  <  00  .  They  are  thus  sequences 
in  a  sequence  space,  say  Y  ,  which  is  a  linear  space.  The  definitions 
of  addition  and  multiplication  by  scalars  in  this  space  are  as  follows: 

OO  00 

Addition:  If  x*  =  {x  }  and  y*  =  {y  }  ,  then 

n  n=o  n  n=o 

00 

z*  =  x*  +  y*  ,  where  z*  -  {z  }  ,  z  =  x  +  y 

y  nn=on  nn 

Multiplication:  ax*  =  aix  }  ={ax}  ,  a  e  cr  . 

r  n  n=o  n  n=o 

If  furthermore,  r*(t)  and  c*(t)  are  bounded,  then  they 
belong  to  the  space  of  bounded  sequences,  say  F  ,  which  is  a  subspace 
of  Y  .  F  is  a  Banach  space  if  we  define  the  norm: 


x*  =  sup 

II  —  co<  p  <00 


X 


* 


X*  €  F 


(2.20) 
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Hereafter  r*(t),  c*(t)  will  be  written  as  r*  ,  c*  and 

their  norms  as  R*,  C*  if  r*,  c*  e  F  . 

In  the  general  situation  when  r*,  c*  e  Y  ,  the  relation 
between  them  may  be  looked  upon  as  a  transformation  in  the  linear  space 
Y  ,  and  hence  can  be  represented  by  the  functional  equation: 

c*  =  B*r*  (2.21) 

where  B*  is  an  operator  which  maps  the  linear  space  Y  into  itself. 
This  operator  is  said  to  be  linear,  if: 


(a) 

B*(xx*  +  x2*)  = 

B*x1*  +  B*x2*  , 

any  x^*> 

x2*  £  Y  . 

(b) 

B*(ax*)  =  aB*x* , 

any  x*  e  Y  , 

and  any 

a  e  . 

B*  is  then  called  a  linear  operator  on  Y  . 

For  the  discrete-data  system  shown  in  Fig.  2.1,  the  input-out¬ 
put  relationship  described  by  the  equation  (2.18)  can  be  easily  verified 
to  be  a  linear  operation  in  Y  .  Hence  we  can  represent  this  relation¬ 
ship  by  the  functional  equation  (2.21),  where  B*  is  a  linear  operator 
defined  by  equation  (2.18).  B*  will  be  called  the  convolution  summation 
operator. 

Now  suppose,  r*,  c*  e  F  ,  then  B*  becomes  a  linear  operator 
in  a  Banach  space,  and  maps  F  -*  F  .  How  this  concept  is  tied  in  with 
the  stability  of  the  linear  system,  and  further  properties  of  B*  ,  if 
the  system  is  stable,  are  discussed  in  the  next  section. 

2 . 4  Stability  of  Linear  Discrete-Data  Systems. 


The  following  definition  of  stability  will  be  adopted  in  this 


’ 
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discussion:  A  discrete-data  system  is  said  to  be  stable  in  the  BIBO 
sense  if  for  any  bounded  input  the  corresponding  sampled  output  is  bounded, 
and  for  every  real  e  >  0  there  exists  a  real  6(e)  >  0  such  that  when¬ 
ever  two  inputs  r^*,  r r*  are  less  than  6(e)  apart,  i.e. 

||  r^*  -  r^*  ||  <  5(e)  ,  the  corresponding  outputs  c^*,  co*  are  less  than 
e  apart,  i.e.  ||  c^*  -  c^*!!  <  e  . 

Referring  now  to  the  linear  mapping  c *  =  B*r*  ,  the  first 
part  of  the  above  definition  implies  that,  if  the  system  is  stable, 
r*,  c*  belong  to  the  space  of  bounded  sequences;  i.e.  they  belong  to 

the  Banach  space  F  ,  and  B*  :  F  F  .  The  second  part  of  the 

definition  further  demands  that  the  mapping  B*  :  F  -*  F  is  continuous. 
Since  B*  is  linear  and  continuous,  B*  is  bounded;  i.e.  a  constant 
M  <  oo  exists  such  that: 

||  B*r*  ||  <  M||r*||  .  (2.22) 

The  smallest  M  for  which  (2.22)  holds  is  called  the  norm  of  the  linear 
operator  B*  ,  denoted  by  ||B*||  . 

According  to  the  above  definition  ||B*||  possesses  the  follow¬ 
ing  properties: 

||  B*r*  ||  <_  ||  B*||||  r*  ||  (2.23) 

for  any  arbitrary  r*  ,  and  there  exists  an  r^*  for  each  e  >  0 
such  that: 

||  B*r  *||  >  (||B*||  -  e)||  r  *||  .  (2.24) 

£  £ 


An  alternate  expression  for  the  norm  of  B*  is: 
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I|b*||  =  Sup  II  B*r*  || 

II  r*||  <L 

22 

This  can  be  shown  to  be  true  as  follows: 

If  ||  r*  ||  <_  1  ,  then: 

||  B*r*  ||  <_  ||  B*  HI  r*||  £  ||  B* 


and  therefore. 


Sup  ||B*r*||  <■  ||  B* 
r*||  <1 


on  the  other  hand  there  exists  an  r^*  for  every  e  >  0  ,  such 


B*re*  II  >  (||  B*  ||  -  e)||  r£* 


If  we  put 


r  * 

ri*  “  irrs- 


then 


=  lTTTfi 


B*r  * 
e 


>  irr*Ti  (lB*H  ■  e)K* 

11  e  11 

=  l|B*||  -  e  . 


Since  ||r^*||  =  1  ,  we  have: 


Sup  ||  B*r*  |)  ||B*r*||  _>  ||B*||  -  e 

r  *||  <1 


(2.25) 


(2.26) 

that : 


and  therefore: 
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Sup  |B*r*||  >  ||  B*  . 
r*||  <1 


(2.27) 


From  (2.26)  and  (2.27)  it  immediately  follows  that: 


B*||  =  Sup  ||  B*r* 
r*||  <1 


(2.28) 


2.5  Computation  of  |B*|  for  a  Given  Linear  Discrete-Data  System. 


For  computational  purpose  it  is  necessary  to  evaluate  an 
expression  for  the  norm  of  B*  in  terms  of  the  system  parameters  for 
any  given  linear  discrete-data  system,  known  to  be  stable.  To  do  this 
we  go  back  to  the  explicit  time  domain  representation  of  the  linear 
discrete-data  system  given  by  equation  (2.18),  viz.: 


c*  -  l 

n=o 


n 


l  h (kT) r (nT-kT) 
k=o 


6(t-nT) 


(2.29) 


From  (2.29): 


00 

n 

=  Sup 

I 

l  h(kT)r (nT-kT) 

6(t-nT) 

— oo<  t  <00 

n=o 

_k=o 

- 

n 

=  Sup 

l 

h(kT)r (nT-kT) 

o<n<°° 

k=o 

00 

<  l  | h(kT) | 

Sup  |r (nT-kT) | 

k=o 

o<n<°° 

00 

«  l  | h (kT) | 

Sup  |r*| 

k=o 

— oo<  t<00 

00 

-  1  | h (kT) | 

Hr*  ||  • 

k=o 


(2.30) 


' 
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Since 


c*  =  B*r* 


B,v  = 


Sup  ||  B*r* 
r*||  <1 


Sup  ||  c*  ||  . 
r*||  <1 


From  (2.30)  letting  ||r*||  =  1  , 


Sup  ||  c*  ||  <  J  |  h(kT) 
r*||  <1  k=o 


i.e. 


B*  ||  1  l  |  h  (kT) 
k=o 


A  classical  argument  can  now  be  used  to  show  that: 


B*  ||  -  l  |  h  (kT)  |  . 

k=o 


To  do  this  choose  an  element  c(nT)  of  the  sequence  c*  , 


n 


c(nT)  =  l  h(kT) r(nT-kT) 
k=o 


n 


=  l  h(nT-kT)r (kT) 
k=o 


For  this  n  ,  choose  r(kT)  =  Sgn  h(nT-kT)  ■  +  1  .  Then 


n  n 

c(nT)  =  l  | h(nT-kT) |  =  l  |h(kT) 
k=o  k=o 


(2.31) 


(2.32) 


n  large. 


Then 
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n 

c* ||  =  ||B*r*||  >_  l  |  h(kT) 

k=o 


Therefore 


n 

||  B*  ||  =  Sup  |jB*r*||  _>  1  |  h  (kT )  |  ,  for  any  n 

||r*||  <1  k=o 


i.e.  ||B*||  >_  l  | h(kT) 

k=o 


From  (2.32)  and  (2.33)  it  follows  that: 


(2.33) 


B* ||  =  l  | h (kT)  |  .  (2.34) 

k=o 


Hereafter  ||B*||  will  be  denoted  by  H*  .  Given  H(z)  a  method  for 
finding  an  upper  bound  for  H*  is  discussed  below. 

It  will  be  assumed  that  H(z)  is  expressible  as  a  ratio  of 
two  polynomials  in  z  ,  and  the  number  of  zeros  of  H(z)  is  at  least 
one  less  than  the  number  of  poles.  This  is  usually  the  case  in  most 
physical  systems,  and  the  impulse  response  of  such  a  system  will  be 
zero  at  t  =  0  ;  in  other  words  h(kT)  =  0  ,  when  k  =  0  ,  so  that 


H*  =  l  |  h  (kT) 
k=l 


Let 


P(z) 

Q(z) 


n 


l 

i-l 


n-i 

z 


n 


I 

i=o 


n-i 

qiZ 


K 


n 

I 

i-l 


Piz 


n-i 


(z  -  z±) 


H(z) 


n 

n 

i-l 


(2.35) 
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where  z^  are  the  n  poles  of  H(z)  .  Here  we  are  assuming  that  H(z) 
has  only  simple  poles.  It  may,  however,  be  mentioned  that  the  results 
can  be  extended  to  multiple  poles  also. 

Applying  equation  (2.19),  for  k  >  0  : 

k-1 

h(kT)  =  Sum  of  the  residues  of  H(z)z 

at  the  poles  of  H(z)  .  (2.36) 

The  residues  for  the  rth  pole  will  be: 


n 


n-i 


Lim 
z  ->  z 


r  L 


(z-z  )P(z)  k_1 

•  z 


Q(z) 


K  l  Plz 

i=l  r  k-1 

-  .  z 

n  r 

II  (z  -  z  . ) 

i-1  r  1 

i^r 


c  z 
r  r 


k-1 


(2.37) 


k-1 

where  c  is  the  coefficient  of  z  and  is  a  constant  for  all  k  . 

r 

Hence 


r 

00 

I 

I  h  (kr 

k-1 

00 

Y 

n 

Y 

L 

k-1 

L 

r=l 

00 

y 

n 

Y 

L 

k-1 

L 

r-1 

c  z 
r  r 


k-1 


k-1 


Since  we  have  assumed  H(z)  is  the  overall  Z-transfer  function  of  a 
stable  discrete-data  system  \z^\  <  1  . 


20 


oo  n 

l  l 


n 


k- 1  r 

c  I  I  z  |  =  l 


k=l  r=l  r  r  r=l 


1  lcrMzr 

k=l 


k-1 


n 


-  I 


i  1  -  z 

r=  1  1  r 


Hence 


n 

H*  1  l 
r=l 


(2.38) 


The  right  hand  side  of  inequality  (2.38)  gives  an  upper  bound  for  H*  . 


} 


CHAPTER  III 


STABILITY  OF  NONLINEAR  DISCRETE-DATA  SYSTEMS 

3. 1  Introduction. 

In  this  chapter  we  will  develop  sufficient  conditions  for  the 
BIBO  stability  of  a  specific  class  of  discrete-data  systems  with  a 
single  continuous  nonlinearity.  Two  such  systems  will  be  considered: 
in  system  A,  the  nonlinearity  will  be  assumed  to  be  situated  in  the  feed¬ 
back  path,  and  in  system  B  it  will  be  assumed  to  be  in  the  forward  path. 
Since  the  method  of  analysis  of  both  systems  is  almost  similar,  the 
analysis  of  system  A  will  be  done  in  detail,  while  that  of  system  B 
will  be  done  more  briefly.  In  either  case  the  boundedness  of  the  out¬ 
put  is  studied  only  at  the  discrete  instants  of  sampling.  The  general 
method  of  approach  is  to  first  express  the  discrete  output  of  the  systems 
as  a  mapping  of  a  Banach  space  into  itself,  and  obtain  sufficient 
conditions  for  the  existence  of  a  fixed  point  within  a  bounded  region 
in  the  space  making  use  of  the  contraction  mapping  theorem  of  Functional 
Analysis.  The  theorem  will  be  stated  and  proved  as  the  analysis  of 
system  A  is  developed. 

3 . 2  System  A:  Description  and  Assumptions. 

The  configuration  of  system  A  is  shown  in  Fig.  3.1.  It  is  a 
single-loop  negative  feedback  system  with  a  single  nonlinear  element  N 
in  the  feedback  path.  The  sampler  is  situated  in  the  forward  path  and 
has  a  constant  period  of  T  .  It  is  followed  by  a  linear  plant  G  , 


, 
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which  may  be  made  of  continuous  controlling  devices,  continuous  and/or 
discrete  compensating  networks  as  well  as  a  zero  order  hold.  r,  c,  and 
e  denote  the  input,  output  and  error  signals  respectively;  they  are 
real  valued  continuous  functions  of  time  t  .  r*,  c*  and  e*  denote 

the  corresponding  discrete  sampled  functions.  r  is  assumed  to  be 
bounded. 


Fig.  3.1  Nonlinear  Discrete-Data  System  A 

The  nonlinear  gain  function  cp ( c)  is  assumed  to  be  continuous 
and  further  satisfies  the  condition  (0 )  =  0  ,  which  permits  its 
representation  in  the  form: 

N  -r 

<J>(c)  -  K  C  +  l  K  C  (3.1) 

i=2 

where  K^,  K_^  are  constants,  and  N  a  finite  integer.  Thus  <Kc)  may 
be  looked  upon  as  consisting  of  a  linear  part  proportional  to  c  ,  and 
a  nonlinear  partcontaining  all  the  higher  order  terms  in  c  .  If  the 
higher  order  terms  were  not  present  the  system  reduces  to  a  linear 
system  with  linear  gain  in  the  feedback  path;  the  system  will 

then  be  said  to  have  been  linearized.  The  BIBO  stability  conditions 
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to  be  derived  for  the  nonlinear  system  will  be  based  on  the  assumption 
that  the  linearized  system  is  BIBO  stable. 

3. 3  Functional  Representation  of  System  A. 

Let  G(s)  be  the  transfer  function  of  the  linear  plant  G  , 
and  G(z)  the  corresponding  Z-transfer  function.  Consider  the  linearized 
system  first.  It  is  shown  in  Fig.  3.2 


r  * 

*  Co 


Fig.  3.2  Linearized  System  Corresponding  to  System  A 

Let  cq  be  the  continuous  output  of  the  linearized  system 
for  the  input  r  ,  and  cq*  be  the  corresponding  sampled  output. 

Also  let  H(z)  be  the  overall  Z-transfer  function  of  the  linearized 
system.  Then  it  can  be  easily  deduced  from  Z-transform  theory  that: 

H(z)  ~  1  TGk'g(z)  (3-2) 

and 

CQ(z)  =  H(z) r (z) 


(3.3) 
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Inverse  transforming  equation  (3.3)  into  the  time  domain  and  using 
functional  notation: 


c  *  =  B*r*  (3.4) 

o 

where  B*  is  the  convolution  summation  operator.  Since  we  have  assumed 
that  the  linearized  system  is  stable,  cq*  is  bounded  whenever  r* 
is  bounded,  and  B*  is  a  bounded  operator;  i.e.  ||B*||  is  finite. 

||B*||  can  be  determined  using  the  formula: 


B* ||  =  H*  =  l  | h (kT)  |  (3.5) 

k=o 


where 


h (kT)  =  Z  ‘Vz) 

H(z)  being  defined  by  equation  (3.2). 


(3.6) 


Also  let 


* 

\ 

-  0  for  t  <  0 

X* 

Sup  |  x*  |  <  °° 

00 

— oo<  r  <°o 

=  \  x(nT)6(t-nT)  for  t  >  0 

.  n=o 

4 

This  is  a  Banach  space  (essentially  )  if: 

||  x*  ||  =  Sup  |x*| 

—  00<  £  <00 


(3.7) 


(3.8) 


In  this  setting,  for  our  stable  linearized  system,  whenever  r*  e  F  , 

c  *  e  F  ;  and  B*  :  F  ■+  F  .  In  other  words  B*  is  a  mapping  of  the 
o 

Banach  space  F  into  itself. 
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Now  consider  the  nonlinear  system  A.  Referring  back  to  Fig.  3.1, 
it  is  clear  that  the  error  signal  e  is  given  by  the  relation: 

N 

e  =  r  -  (c)  =  r  -  K..C  -  J  K.c1  .  (3.9) 

i=2  1 


The  input  to  G  is  the  sampled  error  signal  e*  ,  and 


N 


e*  =  r*  -  4>*(c)  - 


r*  -  K^c* 


-  I 


i=2 


K.c* 

l 


(3.10) 


Laplace  transforming  equation  (3.10): 


e*(s)  =  r*(s)  -  K1c*(s)  - 


N 

l 

1=2 


K±c*J 


(s) 


(3.11) 


Rewriting  equation  (3.11)  in  Z-transform  notation: 


e(z)  ■  r(z)  -  K^c (z)  - 


N 

I 

i=2 


K  c* 


(z) 


(3.12) 


Multiplying  both  sides  of  equation  (3.12)  by  G(z)  ,  and  remembering 
that: 


e(z)G(z)  =  c(z)  and 


G(z) 


1  +  K1G(z) 


=  H  (z) 


the  following  expression  for  c(z)  may  be  derived: 


c ( z )  =  H(z) r (z)  -  H(z) 


N 


i=2 


K.c* 

l 


(z) 


(3.13) 


Transforming  equation  (3.13)  into  the  time  domain  and  using  the  operator 
notation: 


N  ± 

c*  =  B*r*  -  B*  y  K.c* 

. L  n  l 


i=2 


(3.14) 
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with  r*  fixed,  equation  (3.14)  defines  a  relation  of  the  form 


c*  =  Ac* 


(3.15) 


where  A  is  a  nonlinear  operator,  such  that 

v  i 

Ac*  =  B*r*  -  B*  l  K  c*  .  (3.16) 

i=2  1 

3 . 4  Conditions  for  the  BIBO  Stability  of  System  A. 

Whether  the  system  represented  by  equation  (3.14)  is  BIBO  stable 

is  equivalent  to  asking  under  what  conditions  does  it  have  a  solution 

c*  belonging  to  F  ,  if  r*  e  F  ;  or  when  does  A  have  a  fixed 

point  in  F  for  each  r*  e  F.  To  answer  this  question,  we  resort  to 
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a  contraction  mapping  theorem  ,  which  is  stated  and  proved  next.  The 
statement  of  the  theorem  is  slightly  modified  to  suit  our  purpose,  but 
the  proof  is  essentially  the  same  as  that  given  in  reference  23. 


THEOREM:  Let  F  be  a  Banach  space,  and  let  A  be  a  mapping 


of  F  into  F  .  Let  S  be  a  sphere  contained  in  F  of  centre  x^  , 
and  radius  a  ,  such  that: 


(i) 


(3.17) 


for  every  x^,  x^  £  S  ,  and  some  A  :  0  <  X  <  1 


(ii)  || Axq  -  xq||  <  (1  -A  ) a 


(3.18) 


then  there  is  a  unique  fixed  point  x  in  S  ,  that  is  a  point  for 


which  Ax  =  x  . 


Condition  (i)  will  be  called  the  contraction  condition,  and 
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condition  (ii)  the  fixed  point  condition. 


PROOF:  Define  the  sequence  {x^}  which  is 


by: 


x  Ax  -  >  n  1,2,... 
n  n-i 


(and  in  which  xq  is  the  vector  of  our  hypothesis) . 
xq  and  x^  lie  in  the  sphere  S  .  If: 


"  X0II  <  (1  "  *  )a  ,  j  =  1,2,. ..,n 


(so  that  x, ,  x0,  ...,  x  lie  in  S) ,  then: 
1  2*  n 


x  . ,  -  x  ||  <  ||  x  .  ,  -  x  ||  +  ||  x  -  x 


n+1  o 11  —  11  n+1  n 


n 


<  II  Ax  -  Ax  .  II  +  (1  -  An)i 
—  11  n  n-1" 


But 


Ax  -  Ax  ,  <  X  x  -  x  . 

n  n-1 11  —  11  n  n-11 


=  X||Ax  -  -  Ax  _ 
11  n-1  n-2 


,  2.  | 

<  X  x  i  —  x  0 
—  11  n-1  n-2 


<  Xn||x  -  x 
—  11  1  o 


n 


<  X  (1  -  X)a  . 


It  follows  that: 


x  ,  -  -  X  ||  <  xn(l  -  X)a  +  (1  -  Xn)a 
n+1  o"  — 


,n+L 

=  (1  -  X  ) a 


given  recursively 

(3.19) 

By  hypothesis 

(3.20) 

(3.21) 


(3.22) 


(3.23) 


Ir‘:  i 
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so  that  lies  in  S  for  every  n  .  We  see  also  that  because 


ii  .nil 

x  . ,  -  x  <  X  x..  -  x 
n+1  n"  —  11  1  o 


(3.24) 


we  have 


n+r 


n+r-1 

X  II  <  l  ||  x .  ,  -  X. 

n  —  ,L  11  j+1  i 

j=n  J  J 


<  I  x. , ,  -  x . 
—  .  1  J+1  J 

J=n  J  J 


X 


n 


,  x.  -  x 

—  1  -  X  11  1  o 


(3.25) 


For  r  _>  1  ,  the  sequence  {x^}  is  a  Cauchy  sequence,  x^  ->■  x  as 
n  -*  00  ,  and  x  belongs  to  S  . 

Since  Ax  is  continuous  on  S  (by  contraction  condition, 

which  is  a  Lipchitz  condition)  we  now  let  n  -+  00  in  x  =  Ax  -  and 
r  n  n-1 

get 


x  =  Ax 


(3.26) 


There  cannot  be  a  second  point  y  in  S  for  which  y  =  Ay 
for  if  there  were  we  would  have: 

||  x  -  y  ||  =  ||  Ax  -  Ay  ||  <_  X||  x  -  y||  <  ||  x  -  y  ||  (3.27) 

which  is  impossible.  The  proof  is  now  complete. 

In  order  to  apply  the  above  theorem  to  our  nonlinear  mapping 
A  ,  we  must  first  make  sure  that  A  maps  the  Banach  space  F  into  F  . 
This  means  that  in  equation  (3.16)  defining  the  mapping  A  ,  Ac*  should 
belong  to  F  ,  if  r*  and  c*  belong  to  F  .  That  it  is  indeed  so 
can  be  shown  as  follows: 


, 
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Let 


N 

Ac*  =  B*r*  -  B*  l  K 

i=2 

=  B*r*  -  B*Tc* 


(3.28) 


where 


T  :  c* 


N 


i=2 


K..C* 


It  has  already  been  shown  that  B*  :  F  -*  F  ,  if  r*  e  F  . 

Also  T  :  F  F  ,  if  c*  e  F  ,  since 

N  .  N 

!|Tc*||  =  ||  l  K  c*1!  <  l  |K  IHc*!1  <  »  (3.29) 

1=2  1=2 

Therefore  it  follows  that  A  :  F  -*  F  ,  if  r*,  c*  e  F  . 

Next  we  have  to  choose  a  sphere  S  in  F  ,  with  centre  x 

o 

and  radius  ’a'  such  that  the  contraction  and  fixed  point  conditions 
are  satisfied  in  S  .  To  do  this  we  arbitrarily  pick  x^  as  the  out¬ 
put  c^*  of  the  linearized  system  for  an  input  r*  £  F  ;  that  is  we 

let  x  =  c  *  =  B*r*  ,  and  proceed  to  determine  the  conditions  to  be 
oo 

satisfied  by  r*  and  'a'  ,  in  order  that  the  contraction  and  fixed 

point  conditions  are  satisfied  in  S  .  The  choice  of  c  *  as  the  centre 

o 

of  the  sphere  S  is  justified  on  the  basis  that  the  nonlinear  system  A 
may  be  looked  upon  as  a  perturbed  version  of  the  linearized  system,  and 
hence  if  it  has  a  fixed  point  for  a  given  r*  e  F  ,  that  fixed  point 
is  likely  to  be  in  a  neighbourhood  'a*  of  the  solution  of  the  linearized 
system  for  the  same  input  r*  £  F  . 

Consider  now  the  contraction  condition.  Let  c^*,  c2* 


1 
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any  two  points  belonging  to  S  .  Then  the  contraction  condition  may  be 
stated  as: 


Ac1*  -  Ac2*||  <_  X||c1*  -  c2*j|  ,  X  :  0  <  X  <  1 


(3.30) 


From  equation  (3.16): 


Ac1*  -  Ac2* 


N 


=  ||  B*  l  K.c,,*1  - 


N 


i=2 


B*  l  K.c 
1=2 


N 


ills*  II  l  IlK.Cc^1 

i=2 


-  c. 


-  H*  l  ||K.(c2*  -  Cj*)  l  c1*1  1  k  c2*k| 
i=2  k=o 


N 


i-1 


<  H*||c2*  -  Cl*||  l  |K±|  l  ||c  * 

i=2  k=o 


i-l-kn  *„k 
c2  11 


(3.31) 


Since  c^*,  c2*  €  s  : 


\\c±*\  >  Hc2*l  1  II  C05V II  +  a 

=  ||  B*r*  ||  +  a 
1  ||  B*  1  Hr*  ||  +  a 

=  H*R*  +  a  (3.32) 


II  aC;l*  -  Ac2*|| 

N  i-1  ,  ,  , 

£H*||c2*  -  c1*||  l  | K.  |  l  (H*R*+a)  (H*R*+a) 

i=2  k=o 

N 

£H*||c2*  -  Cl*\\  l  1^1  i(H*R*+a)1  1  .  (3.33) 
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Hence  (3.30)  will  be  satisfied  if: 


N 


H*||  c2* 


-  c 


*||  l  |K.  |i(H*R*+a)X  1  A||c 
i=2  1 


*  -  r  * 


Since  ||  c2*  -  c^*||  =  ||  c^*  -  c2*||  ,  (3.34)  is  equivalent  to: 


N 


H*  l  |K.  |  i(H*R*+a)1  1  £  X  . 
i=2  1 


Now  consider  the  fixed  point  condition.  With  c  *  as 

o 

centre  of  the  sphere  S  ,  this  may  be  stated  as  follows: 


Ac  *  -  c  *  <  (1  -  A)  a  . 

o  o  11  — 


Since 


c  *  =  B*r* 
o 


c  *  <  H*R*  . 

o  "  — 


Also  it  follows  from  equation  (3.16)  that: 


N 


Ac  * 
o 


=  c 


*  _ 


o 


B*  y  K.c  *x 

1  O 

1=2 


Transposing  cq*  to  the  left  hand  side  and  taking  norms  on  both 


N 


Ac  *  -  c  *  = 

o  o  11 


-B*  y  K.c  *x 

i=2  1  0 


N 


<  ||b*||  y  Ik.  I  lie  *i 

-  »  ia2  1  i1"  o 


N 


_<  H*  l  |K.  |  (H*R*)  . 

i=2  1 


(3.34) 

(3.35) 

the 

(3.36) 

(3.37) 

(3.38) 

(3.39) 

sides : 


(3.40) 


. 
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Hence  (3.36)  will  be  satisfied  if: 

N 

H*  l  | K. |  (H*R*)X  <  (1  -  X)a  .  (3.41) 

i=2  1 

Recapitulating,  given  an  aribtrary  input  r*  e  F  ,  with 
||  r* ||  =  R*  ,  inequalities  (3.35)  and  (3.41)  give  the  sufficient  conditions 
for  the  nonlinear  mapping  expressed  by  equation  (3.16)  to  have  a  fixed 
point  c*  in  a  neighbourhood  'a'  of  the  solution  of  the  stable 
linearized  system  for  the  same  input  r*  .  Moreover,  since  c*  belongs 
to  the  sphere  S  ,  an  upper  bound  for  c*  is  given  by: 

||  c* ||  1  H*R*  +  a  .  (3.42) 

We  are  now  in  a  position  to  state  the  following  theorem  which 
embodies  sufficient  conditions  for  the  existence  of  a  region  in  which 
the  system  A  exhibits  bounded- input  bounded-output  stability: 

THEOREM  1:  The  system  A,  satisfying  the  conditions  set  forth 
in  section  3.2,  and  driven  by  an  arbitrary  input  r  ,  with  ||r*||  =  R*  , 
possesses  a  bounded  sampled  output  c*  in  a  neighbourhood  ’a'  of  the 
solution  of  the  stable  linearized  system  for  the  same  input  r  ,  if 
a  X  :  0  <  X  <  1  ,  exists  such  that  the  following  inequalities  are 
satisfied : 

N 

H*  l  | K . | i(H*R*+a)X  ^X  (3.43) 

i=2  1 


H*  l  | K. | (H*R*) 1  <  (1  -  X)a  . 
i=2  1 


(3.44) 


33 


An  upper  bound  for  c*  is  then  given  by: 

||  c* ||  £  H*R*  +  a  .  (3.45) 

3. 5  Conditions  for  the  BIBO  Stability  of  System  B. 

System  B  differs  from  system  A  in  that  the  nonlinearity  is 
situated  in  the  forward  path  instead  of  in  the  feedback  path.  The 
configuration  of  system  B  is  shown  in  Fig.  3.3,  and  that  of  the 
corresponding  linearized  system  in  Fig.  3.4. 


/  \ 


Fig.  3.3  Nonlinear  Discrete-Data  System  B 


Fig.  3.4  Linearized  System  Corresponding  to  System  B 


The  equation  for  the  gain  is: 


N  . 
4> (e)  -  K^e  +  l  K^e 


i-2 


(3.46) 
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where 

e  =  r  -  c  (3.47) 

The  linearized  system  is  assumed  to  be  stable  as  before.  The  overall 
Z-transfer  function  of  the  linearized  system  is: 


H(z)  = 


K^Cz) 


1  +  KjGCz)  * 


(3.48) 


The  equation  for  the  nonlinear  system  can  be  derived  as  follows: 


c (z )  =  G( z)<p(e)  (z) 


N 


-  G (z)  [K  e(z)  +  (  l  K  e*1) (z) ] 

i=2 


N 


=  K  G(z)[r(z)  -  c  ( z)  ]  +  G  (z)  (  £  K  e*1)  (z)  .  (3.49) 

i=2 


Transposing  -K^G(z)c(z)  to  the  left  hand  side,  and  dividing  both  sides 


by  1  +  KxG(z)  : 


K  G(z)r(z)  K  G(z)  N 

=  m^(7)  +  *  1  +  KxG(z)  (  Ki6*  )(z) 


N 


H(z)r(z)  +  ~~  H (z) (  l  K  e*i)(z)  . 

K1  i=2  1 


(3.50) 


Inverse  transforming  (3.50)  and  using  operator  notation: 


1  v  i 
c*  =  B*r*  +  -f—  B*  l  Ke* 
K,  .  i 


1  i=2 


(3.51) 


where 


e*  =  r*  -  c* 


(3.52) 
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With  r*  fixed,  equation  (3.51)  defines  a  relation  of  the  form: 

c*  =  Ac*  (3.53) 

where  A  is  a  nonlinear  operator,  such  that: 

N 

Ac*  =  B*r*  +  ■—  B*  y  K.e*1  .  (3.54) 

K1  i=2  1 

It  is  easy  to  show  as  before  that  A  maps  a  Banach  space  F 
into  F  if  r*,  c*  £  F  . 

To  investigate  the  existence  of  a  fixed  point  for  the  mapping 

A  in  a  sphere  S  with  centre  c  *  =  B*r*  ,  and  radius  ’a'  consider 

o 

first  the  contraction  condition.  Let  c^*,  c2*  be  anY  two  P°ints  in 
the  sphere  S  .  Let  e^*,  e2*  be  t*ie  sampled  error  functions 
corresponding  to  c^*,  c2*  * 

From  equation  (3.54): 


||  Ac1*  -  Ac2* 


N  .  N 

I  K.e^1  -  B*  I  k^*1 

1=2  1=2 


N 


I  ||K  (e  A1  -  e  *5 
1=2  1  1 


N 


=  l 


1=2 


W 


i-i 


I 

k=o 


i-l-k 


N  i-1 


II  I  |KJ  l 

1=2  k=o 


i-l-k 


(3.55) 


Note:  Since  the  linearized  system  is  assumed  to  be  stable,  has  got 

to  be  positive;  hence,  the  absolute  value  sign  is  not  put  on  . 


} 
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Now  for  a  given  r*  it  follows  from  equation  (3.52)  that: 


e  *  —  e  *  =  r  *  —  r  * 

1  e2  11  11  1  c2 


(3.56) 


el*  ||  £  R*  +  II  c1* 


(3.57) 


e2*||  1  R*  +  ||  c2* 


(3.58) 


Since  c^*,  c2*  €  S  ,  we  also  have 


cx*  || ,  f|  c2*  (I  <_  H*R*  +  a  . 


(3.59) 


Therefore 


e1*  ||  <_  R*  +  H*R*  +  a 


e^*  ||  R*  +  H*R*  +  a 


(3.60) 

(3.61) 


Making  use  of  (3.56),  (3.60)  and  (3.61)  we  may  write: 


Ac^*  -  Ac2* 


N 


i=2 

N 


1  «  i-1  ,  ,  ,  , 

<_  £-H*||Cl*  -  c2*||  l  |K  I  l  (R*+H*R*+a) 1  (R*+H*R*+a) 

'  ‘  "  k=o 


H* 


K  ~  „c1*  -  c2*||  l  |K.|  i(R*+H*R*+a) 
1  i=2 


i-1 


(3.62) 


In  order  to  satisfy  the  contraction  condition  we  require: 


Ac^*  -  Ac2*||  <_  X||  c1*  -  c2*  ||  ,  X  :  0  <  X  <  1  . 


(3.63) 


By  virtue  of  (3.62),  (3.63)  will  be  satisfied  if: 


1  N  i-1 

~H*||Cl*-  c2*||  l  |Ki|i(R*+H*R*+a)X  £  X||  C;L*  -  c2* 


i=2 


(3.64) 


? 
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or  if : 


1  N  -_i 

77—  H*  l  |K.  | i (R*+H*R*+a) 1  <_  A  . 

1  i=2  1 


(3.65) 


Now  consider  the  fixed  point  condition,  which  requires  that: 


Ac  *  -  c  *  <  (1  -  A)a  . 

o  o  11 


(3.66) 


If  eQ*  is  the  sampled  error  corresponding  to  cq*  ,  it 
follows  from  equation  (3.54)  that: 


N 


Ac  * 
o 


c  *  +  ~~  B*  y  K.e  *i 
0  K1  i=2  1  ° 


(3.67) 


Transposing  cq*  to  the  left  hand  side  and  taking  norms  on  both  sides: 


Ac  *  -  c  * 
o  o 


1  N  * 
~  «k7  b*  l  KiV 


'1  i=2 


i  3  |  in  Ni 

<  “  H*  J  K.  e  *  . 

—  K.  .  1  l 111  o  11 


1  i=2 


(3.68) 


But 


e  *  <  R*  +  c  * 

o  —  11  o  11 

<  R*  +  H*R*  . 


(3.69) 


Therefore 


Ac  *  -  c  *||  £  ~  H*  l  |K.  |  (R*  +  H*R*) 1  . 

0  °  ±=2  1 


N 


(3.70) 


Hence  the  fixed  point  condition  will  be  satisfied  if: 


i  N  i 

—  H*  l  |K.|(R*  +  H*R*)  <  (1  -  A)  a  . 

K.  ,  o  i 


1  i=2 


(3.71) 
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Thus  inequalities  (3.65)  and  (3.71)  are  sufficient  for  the 
mapping  A  to  have  a  fixed  point  c*  inside  the  sphere  S  .  Since 
c*  lies  inside  S  ,  an  upper  bound  for  c*  is  given  by: 

||  c*||  £  H*R*  +  a  .  (3.72) 

We  may  now  state  the  following  theorem  concerning  the  bounded- 
input  bounded-output  stability  of  system  B. 

THEOREM  2:  The  system  B,  conforming  to  the  description  and 
assumptions  set  forth  at  the  beginning  of  this  section,  and  driven  by 

an  arbitrary  input  r  ,  with  |jr*j|  =  R*  ,  possesses  a  bounded  sampled 

output  c*  in  a  neighbourhood  'a*  of  the  solution  of  the  stable 
linearized  system  for  the  same  input  r  ,  if  a  X  :  0  <  A  <  1  ,  exists 
such  that  the  following  inequalities  are  satisfied: 

1  N  '-1 

——  H*  y  |K. I i(R*  +  H*R*  +  a)1  <  X  (3.73) 

K-i  .  o  i 

1  i=2 

'  N 

H*  l  | K . | (R*  +  H^R*)1  <  (1  -X)a  .  (3.74) 

K1  i=2  1 

An  upper  bound  for  c*  is  then  given  by: 

||  c* II  <  H*R*  +  a  .  (3.75) 


CHAPTER  IV 


COMPUTATION  OF  INPUT  AND  OUTPUT  BOUNDS 

4.1  Introduction. 

In  the  last  chapter  we  formulated  the  sufficient  conditions 
for  the  existence  of  a  bounded  output  for  two  different  types  of  non¬ 
linear  discrete-data  systems.  In  practice  we  are  interested  in 
determining  the  maximum  value  of  R*  for  which  the  sytems  exhibit  a 
bounded  output.  To  do  this  we  have  to  choose  the  optimum  values  of 
’a’  and  X  which  give  the  maximum  value  of  R*  satisfying  the 
contraction  and  fixed  point  conditions.  This  has  to  be  done  only  by 

trial  and  error  in  the  general  case  of  the  nonlinearity.  However,  if 

th 

the  nonlinear  part  of  the  gain  function  consists  of  a  single  n 
degree  term  a  direct  approach  is  available  to  determine  the  maximum 
value  of  R*  to  ensure  a  bounded  output,  and  the  corresponding  upper 
bound  for  the  output.  This  method  will  be  described  in  the  following 
sections,  together  with  numerical  examples  to  illustrate  the  computational 
procedure . 

th 

4.2  Input  and  Output  Bounds  for  System  A  with  a  Single  n  Degree 

Term  in  the  Nonlinearity. 

The  above  system  differs  from  the  general  system  in  that  the 
nonlinear  gain  function  is  given  by: 

<j>(c)  =  K  c  +  K  c”  ,  n  >_  2  .  (4.1) 

The  equation  for  c*  will  then  be: 


f '  «! 
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c*  =  B*r*  -  B*K  c* 

n 


n 


The  contraction  condition  reduces  to: 


nH* I K  I (H*R*  +  a)n  1  <  X 
i  n  1  — 


and  the  fixed  point  condition  to: 


H*|K  I (H*R*)n  <  (1  -  X) a  . 
n 


Note  that  X  :  0  <  X  <  1  ,  and  a  >  0  .  Letting 


a  =  bH*R* 


inequality  (4.3)  may  be  written  as: 


nH*|K  | (H*R*  +  bH*R*)n  1  £  X 


or 


nH* I K  I (H*R*)n_1(l  +  b)n  1  <  X 
n  — 


or 


Ik  I  H*(H*R*)n  1  <  - - - - 

n  /  -1  i  i_  \  ^ 

n(l  +  b) 


Similarly  inequality  (4.4)  may  be  written  as: 


H*|K  |(H*R*)n  <  (1  -  X)bH*R* 
n 


or 


Let 


Then  (4.6)  implies 


K  I H* (H*R*)n  1  <  (1  -  X)b  . 
n 


| K  | H* (H* (H*R*) n  1  =  x  . 


x  < 


X 


n(l  +  b) 


n-1 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


(4.7) 

(4.8) 

(4.9) 
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(4.7)  implies 


x  <  (1  -  A)b  .  (4.10) 

To  find  the  maximum  value  of  x  ,  satisfying  (4.9)  and  (4.10) 
note  that  for  a  given  n  ,  and  b  >  0  ,  the  right  hand  side  of  (4.9) 
is  an  increasing  function  of  A  ,  whereas  the  right  hand  side  of  (4.10) 
is  a  decreasing  function  of  A  .  Hence  for  a  given  n  and  b  ,  x  will 
be  a  maximum  for  that  value  of  A  at  which  the  two  curves: 


y  = 


n(l  +  b) 


n-1 


and 


(1  -  A)b  intersect  . 


Similarly  for  a  given  n  and  A  ,  0  <  A  <  1  ,  the  right  hand 

side  of  (4.9)  is  a  decreasing  function  of  b  ,  whereas  the  right  hand 
side  of  (4.10)  is  an  increasing  function  of  b  .  Hence  for  a  given  n 
and  A  ,  x  will  be  a  maximum  for  that  value  of  b  at  which  the  two 
curves : 

y  =  -  ■  and  y  =  (1  -  A)b  intersect  . 

n(l  +  b) 

Therefore  when  A  and  b  are  both  varied  for  a  given  n  ,  the 
maximum  value  of  x  will  lie  on  the  curve  of  intersection  between  the 
two  surfaces: 

y  =  - X  and  y  -  (1  -  X)b  , 

n(l  +  b)"'1 

and  will  correspond  to  the  maximum  height  of  that  curve,  from  the 
A-b  plane. 

To  obtain  y  proceed  as  follows: 
max 
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At  the  point  of  intersection  of  the  two  surfaces: 


X 

n(l  +  b)n_1 


(1  -  A)b 


from  which: 


nb (1  +  b)n  1 

nb (1  +  b)n_1  +  1 


(4.11) 


Note  that  from  (4.11)  it  is  obvious  that  the  condition  0  <  X  <  1  is 
satisfied  for  a  given  n  and  b  >  0  .  Also 


y  =  (1  -  A)b 

_ b _ 

nb (1  +  b)n_1  +  1 


(4.12) 


The  condition  for  y  to  be  maximum  are: 

*  *  0  ’  (4-13) 

and 

,2 

^  <  0  ,  (4.14) 

db 

for  the  value  of  b  for  which  (4.13)  holds. 

Differentiating  (4.12)  with  respect  to  b  ,  and  equating  to 
zero  we  get: 


1  -  n(n  -  l)b2(l  +  b)n  2  =  0 


(4.15) 


The  positive  real  root  of  equation  (4.15)  gives  the  value  of 
b  for  which  y  is  a  maximum.  It  may  be  observed  that  the  left  hand 
side  of  equation  (4.15)  is  a  decreasing  function  of  b  ,  for  b  >  0  , 
and  since  it  is  equal  to  1  when  b  =  0  ,  the  equation  possesses  one 


1 

i  oil  (t  i  i)  b  i 
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and  only  one  positive  real  root.  Moreover,  it  can  be  easily  verified 
that  for  this  root  the  second  derivative  of  y  with  respect  to  b  is 
less  than  zero.  Substituting  this  root  in  equation  (4.12)  we  get  the 
maximum  value  of  y  . 

Since  Ymax  represents  the  maximum  value  of  x  ,  satisfying 
the  inequalities  (4.9)  and  (4.10),  and  x  is  given  by  equation  (4.8) 
we  may  conclude  that  the  maximum  value  of  R*  satisfying  both  the 
contraction  and  fixed  point  conditions  is  given  by  the  inequality: 


K  jH*(H*R*)n~1  <  - 2 — r - 

n  nb  (1  +  b  )n  1  +  1 

o  o 


(4.16) 


where  bQ  represents  the  positive  real  root  of  equation  (4.15).  From 
(4.16)  it  follows  that  the  best  possible  bound  on  R*  for  the  system 
to  exhibit  a  bounded  output  is  given  by  the  following  inequality: 


R*  < 


H* 


IK  !H*{nb  (1  +  b  )n"1  +  1} 
'  n 1  o  o 


1 

n-1 


(4.17) 


The  corresponding  bound  on  c*  is  then  given  by: 


C*  <  H*R*  +  a 


=  H*R*(1  +  b) 


<  (i  +  b  > 


b 

o 


K  |H*{nb  (1  +  b  )n"1  +  1} 
n 1  o  o 


1 


n-1 


(4.18) 


Numerical  Example  1: 


As  a  numerical  example  consider  the  nonlinear  discrete-data 


(21 
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system  shown  in  Fig.  4.1. 


Nonlinearity 
<t>(c)  =  c+K^c3 


Fig.  4.1  Discrete-Data  System  with  a  Cubic  Nonlinearity 

in  the  Feedback  Path. 


It  has  a  cubic  type  of  nonlinearity  represented  by  the  gain  function: 

<j>(c)  =  c  +  K3c3  .  (4.19) 

Referring  to  equation  (4.1)  this  corresponds  to  the  case  when  =  1  , 

and  n  =  3  .  The  forward  path  of  the  system  contains  a  sampler  with 

period  T  =  1  second,  followed  by  a  linear  plant  consisting  of  a  zero 

order  hold  and  a  controller  with  transfer  function  — y-r  ,  so  that  the 

s+1 

Z-transfer  function  of  the  linear  plant  may  be  written  as: 


G(z) 


_ 1 

s  + 


0.632 
z  -  0.368 


(4.20) 


and 
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H(z)  = 


G(z) 


1  +  G  (z) 

0.632 
z  +  0 . 264 


(4.21) 


Inverse  transforming  H(z)  : 

h(kT)  =  £  Residues  of  H(z)z 


^  ^  at  the  poles  of  H(z)z^  ^  .  (4.22) 


Inserting  numerical  values  and  evaluating  the  residues 


h  (kT) 


Hence 


for  k  =  0 

.632  (- 

0.264)^  for  k  >  0 

00 

=  V 

L 

1  h (kT)  | 

k=o 

OO 

-  1 

| 0.632  (-0 . 264)k_1| 

k=l 

«  0.632 

OO 

l  0.264 

k=l 

0.86 


Now  for  n  =  3  ,  equation  (4.15)  reduces  to: 


1  -  6b  (1  +  b)  =  0 


(4.23) 


(4.24) 


(4.25) 


(4.25)  is  a  cubic  equation,  solving  which  we  find  that  it  has  a  positive 
real  root  =  0.35  .  Hence  bQ  =  0.35  .  Now  letting  H*  =  0.86  ,  n  =  3 
and  bQ  =  0.35  in  inequality  (4.17)  the  maximum  permissible  sampled 
input  to  ensure  the  existence  of  a  bounded  sampled  output  will  be 


given  by: 
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R*  < 


0.86 


0.35 


1 

2 


=  0.435  K 


i |  0 . 86 {3  x  0.35  x  1.352  +  l}_ 

1 

2 


(4.26) 


To  obtain  the  corresponding  bound  on  the  sampled  output  use 
inequality  (4.18)  and  obtain: 


C*  <  0.795 


(4.27) 


Note  that  the  optimum  value  of  X  in  this  case  is  0.66  from 

equation  (4.11)  by  letting  n  =  3  ,  and  bQ  =  0.35  . 

Inequality  (4.26)  can  also  be  used  to  find  the  permissible 

maximum  value  of  |K^j  for  a  specified  bcund  on  R*  .  For  example  in 

the  case  of  a  unit  step  input  R*  =  1  ,  and  I  K_  j  works  out  to 

3  1  max 

0,188  ,  and  the  corresponding  bound  on  C*  becomes  1,83  from  (4.27). 

An  example  of  the  same  system,  where  the  solution  obtained  by 

solving  the  difference  equation  of  the  system  by  means  of  a  recurrence 

relation  for  the  case  =  0.1  ,  and  a  unit  step  input,  may  be  found 

in  reference  24.  In  that  example,  c*  occurs  at  the  end  of  the 

max 

first  sampling  period,  and  is  equal  to  0.632  .  Using  the  contraction 
technique  outlined  above  the  bounds  for  R*  and  C*  work  out  to  1.37 
and  2.5  when  =  0.1  ,  which  are  consistent  with  the  results 

obtained  in  reference  24. 
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4 . 3  Input  and  Output  Bounds  for  System  B  with  a  Single  n 

Term  in  the  Nonlinearity. 


th 


For  the  above  system  the  nonlinear  gain  function  is 


<J>(e)  =  K.e  +  Ken,  n>2,  K  >  0 
in  —  i 


and  the  equation  for  c*  is: 


c*  =  B*r*  +  ~  B*  K  e*n 

n 


The  contraction  condition  reduces  to: 


nH* I K  I (R*  +  H*R*  +  a)n  1  <  X 
1  n1  — 


and  the  fixed  point  condition  to: 


—  H*|K  j (R*  +  H*R*)n  <  (1  -  A) a 
n 


Now  we  let 


a  =  b (R*  +  H*R*) 


Then  inequality  (4.30)  may  be  written  as: 


nH* i K  | (R*  +  H*R*)n  1(1  +  b)n  1  <  A 
'  n  — 

or 


K 


—  H* (R*  +  H*R*)n  1  < 


K 


X 


n(l  +  b) 


n-1 


Similarly  inequality  (4,31)  may  be  written  as: 


K 


—  H* (R*  +  H*R*)n  1  <  (1  -  X)b 


K 


Degree 

given  by: 

(4.28) 

(4.29) 

(4.30) 

(4.31) 


(4.32) 

(4.33) 
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Letting 


K 

n 


H*(R*  +  H*R*)n  1 


x 


(4.34) 


(4.32)  implies 


X 

x  <  - T 

n(l  +  b)n  1 

(4.33)  implies 


(4.35) 


x  <  (1  -  X)b 


(4.36) 


Note  that  the  inequalities  (4.35)  and  (4.36)  are  the  same  as 
inequalities  (4.9)  and  (4.10)  for  system  A  considered  in  section  4.2, 
though  x  in  this  case  represents  a  different  quantity.  The  maximum 
value  of  x  for  which  the  inequalities  (4.35)  and  (4.36)  are  satisfied 
can  therefore  be  found  as  before,  and  the  condition  for  BIBO  stability 
of  the  system  reduces  to: 

|K  |  b 

_JL_  H*(R*  +  H*R*)n  <  - 2 — — -  (4.37) 

K1  nb  (1  +  b  )n  1  +  1 

o  o 


where  b  is  the  real  positive  root  of  equation  (4.15). 
o 

From  (4.37)  it  can  be  easily  deduced  that  the  best  bound  for 


R*  for  the  existence  of  a  bounded  sampled  output  is  given  by: 

1 


R*  < 


1  +  H* 


K.  b 
1  o 


|k|  H*{nb  (1  +  b  )n  1  +  1} 
L-  1  n  o  o 


n-1 


(4.38) 


The  corresponding  bound  on  C*  is  given  by: 
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C*  ^  H*R*  +  a 

=  H*R*  +  b(R*  +  H*R*) 

1 

n-1 

(4.39) 


<  H*  +  b (1  +  H*)  R* 


H*  ,  , 

1  +  H*  bo 


K1  b 
1  o 


L|K  | H*{nb  (1  +  b  )n  1  +  1} 
'  n  o  o 


Numerical  Example  2: 


As  a  numerical  example  we  will  consider  the  nonlinear 


V 


* 


Fig.  4.2  Discrete-Data  System  with  a  Cubic  Nonlinearity 

in  the  Forward  Path. 

discrete-data  system  shown  in  Fig.  4.2,  It  will  be  observed  that  the 
linear  part  of  the  system  is  the  same  as  that  of  the  system  considered 
in  example  1;  the  nonlinearity  is  however  situated  in  the  forward  path 
instead  of  in  the  feedback  path.  The  nonlinear  gain  function  is  given 
by: 

4>(e)  =  e  +  K3e3  .  (4.40) 

The  period  T  of  the  sampler  is  again  assumed  to  be  1  second.  Since 
the  linearized  system  corresponding  to  the  above  system  is  the  same  as 
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that  of  the  one  considered  in  example  1,  H*  =  0.86  ,  for  this  system 
as  well.  Also  n  =  3  ,  K^  =  1  ,  and  bQ  =0.35  .  Substituting  these 
values  in  inequality  (4.38)  the  maximum  permissible  sampled  input  to 
ensure  the  existence  of  a  bounded  sampled  output  will  be  given  by: 


R*  < 


1 


1.86 


0.35 


JL 

2 


,  K3  |  0 . 86  {3  x  0.35  x  1„  35  +  1 }  J 


0.2  IK 


1 

2 


(4.41) 


The  corresponding  bound  on  the  sampled  output  can  be  obtained 
using  inequality  (4.39).  It  is  given  by: 


c* <  [itH  +  °'35-' 


0.35 


K  1 0.86  (3  x  0.35  x  1.35  +  l}_ 


1 

2 


<  0.478  |K, 


1 

2 


(4.42) 


As  before,  inequality  (4.41;  can  be  used  to  find  the  permissible 

maximum  value  of  |K^|  for  a  specified  bound  on  R*  .  Thus  in  the  case 

of  a  unit  step  input  | K  |  works  out  to  0.04  ,  and  the  correspond- 

3  max 

ing  bound  on  C*  becomes  2.39  from  (4.42). 


, 
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CHAPTER  V 


STABILITY  OF  A  DISCRETE-DATA  SYSTEM  WITH  A 
SLOPE  RESTRICTED  NONLINEARITY 


5 . 1  Introduction. 

In  this  chapter  we  will  present  and  prove  a  criterion  for 
the  bounded-input  bounded-output  stability  of  a  class  of  discrete-data 
systems  with  a  slope  restricted  nonlinearity.  The  proof  is  based  on 
the  principle  of  contraction  mapping.  A  numerical  example  is  worked 
out  to  illustrate  the  application  of  the  criterion,  and  to  discuss  its 
implication  in  relation  to  some  existing  results. 

5 . 2  System  Description. 

The  configuration  of  the  system  considered  is  shown  in 
Fig .  5 . 1(a) . 

a 

Fig.  5.1(a) 

K„e 

e) 


d|(el  c 

de  2 

Fig.  5.1(b) 


V 


oV 


Fig.  5.1  Class  N  Nonlinear  Discrete-Data  System 
(a)  System  Configuration  (b)  Gain  Function. 
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The  nonlinear  gain  function  cf (e )  (shown  in  Fig,  5, 1(b))  is  assumed  to 
be  continuous  and  differentiable,  and  further  satisfies  the  condtions: 

d(0)  =  0  (5,1) 

o  <  K1  e  K2  „  (5.2) 

Note  that  (5.1)  and  (5.2)  imply  that  cp(e)  is  confined  to  the  gain 
sector  [K^,  K^]  i.e.  0  <  K1  <  »  The  rest  of  the  system 

comprises  a  sampler  with  period  T  ,  followed  by  a  linear  plant  with 
transfer  function  G(s)  .  Negative  unity  feedback  is  assumed.  The 
system  so  described  will  be  designated  as  belonging  to  the  class  N  . 

5 .3  Stability  Criterion  for  N  . 

Let  b(e)  be  replaced  by  a  linear  gain  K  ,  and  the 
resulting  linear  system  designated  by  L  ,  Define 

OC 

H*  =  l  |h(kT) |  , 

k=o 

where  h(kT)  is  the  closed-loop  impulse  response  sequence  of  L  . 

Then  the  stability  criterion  for  N  may  be  formulated  as  follows: 

THEOREM  1;  A  nonlinear  discrete-data  system  belonging  to 

the  class  N  is  BIBO  stable,  if  a  positive  K  exists  for  which 

K  K 

L  is  stable,  and  K  -  and  <  K  +  —  . 

PROOF :  Choose  a  K  for  which  L  is  stable;  this  implies 

that,  if  the  input  -  output  relation  of  L  is  represented  by  the  linear 


functional  equation: 


f 
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c*  =  B*r* 


(5.3) 


where  B*  is  the  convolution  summation  operator,  then: 


||b*||  =  h*  < 


(5.4) 


Note  that  in  (5.3)  r*,  c*  denote  the  sampled  input  and  output  respectively, 

and  since  L  is  stable  r*»  e*  belong  to  a  Banach  space  F  ,  equipped 

with  the  uniform  norm  ||x*j|  =  Sup  jx*|  ,  and  (5,4)  defines  the  norm 

x*  c  F 

of  B*  ,  induced  by  the  norm  defined  in  F  . 

Now  consider  the  nonlinear  system  N  .  From  Z-transform  theory: 

c ( z )  =  4> (e)  (z)G(z) 

=  Ke(z)G(z)  +  4> (e )  (z)G(z)  -  Ke(z)G(z) 

=  K[r(z)  -  c(z)  ]G(z)  +  <p(e)  (z)G(z)  -  Ke(z)G(z)  .  (5.5) 

Transposing  -  Kc(z)G(z)  to  the  left  hand  side,  dividing  both  sides  by 
1  +  KG(z)  ,  and  putting  / — r  =  H(z)  ,  (5,5)  reduces  to: 

L  t  Js.Lt  {  Z ) 


c ( z )  *  H(z) r(z)  +  H(z)  _  e(z)J 


(5.6) 


Inverse  transforming  (5.6)  and  using  the  operator  notation: 


c*  =  B*r*  +  B*[—  -  e*] 

K 


(5.7) 


Since  for  a  given  r*  , 


q  /k  m  ^  k  —  £  k 


(5.8) 


(5.7)  defines  a  relation  of  the  form: 
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c*  =  Ac* 


(5.9) 


where  A  is  a  nonlinear  operator,  such  that: 

Ac*  =  B*r*  +  B*[^f -  e*]  .  (5.10) 

K 

Now  if  r*,  c*  e  F  ,  A  maps  F  -+  F  .  This  can  be  shown  as 

follows : 

||  Ac*  ||  =  ||  B*r*  +  B*e*  [j  •  -  1]  || 

1  II  B*  ||  |  r*  ||  +  ||B*||||e*||[|  +  1]  .  (5.11) 

Note  that  in  (5.11): 


(1) 

|  B*  =  H*  < 

00  ,  since  L 

is 

stable 

(2) 

|l  r*  ||  <  -  , 

as  r*  £  F  . 

(3) 

||e*  ||  <  ||  r* || 

+  ||  c*  ||  ,  from 

(5 

.8) 

<  00  ,  as  r*,  c*  £  F  . 


(4)  ||  y— — -|1  <  <  00  ,  which  is  a  consequence  of  the 

hypotheses  (5,1)  and  (5.2). 

Thus  all  quantities  on  the  right  hand  side  of  (5.11)  are 
bounded,  and  hence  jj Ac* j|  <  00  ,  which  means  A  :  F  F  . 

We  can  therefore  apply  the  contraction  mapping  theorem  to 
determine  the  conditions  under  which  the  mapping  A  has  a  fixed  point 
in  F  .  To  do  this  we  choose  the  output  c  *  of  the  linear  system  L 
to  a  bounded  input  r  ,  as  a  first  approximation  for  the  solution 
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of  the  nonlinear  system  N  for  the  same  input  r  ,  and 

show  that  a  sphere  S  with  a  finite  radius  'a'  ,  centred  on  c  * 

o 

exists  in  which  the  fixed  point  condition  is  always  satisfied.  Next 
we  apply  the  contraction  condition  to  the  points  inside  this  sphere, 
and  obtain  the  requirement  for  the  existence  of  a  fixed  point  in  S  . 

Consider  first  the  fixed  point  condition  which  states: 


Ac  *  -  c  *  <  (1  -  A)a  ,  A  :  0  <  A  <  1  .  (5.12) 

o  o 


Since  L  is  stable  c  *  is  bounded;  so  also  is  e  *  the  corresponding 

o  o 

sampled  error.  Also: 


c  *  =  B*r* 
o 


(5.13) 


Hence  from  equation  (5.10) 


Ac  *  -  c  *  =  B* 
o  o 


4>(eo*) 

K 


-  e 


* 


B*e  * 
o 


<f>(eQ*) 

Ke  * 
o 


-  1 


(5.14) 


Therefore 


Ac  *  -  c  *  < 

o  o  "  — 


B*  e 


* 


1 

K 


4>(eQ*) 

e  * 
o 


+  1 


(5.15) 


Note  that  in  (5.15): 


(1)  B*  |  =  H*  <  00  ,  since  L  is  stable, 


(2)  ||e  *||  <  00  ,  since  e^*  is  bounded, 


<Keo*) 


<  k2  < 


(3) 


00 


,  which  is  a  consequence  of  the 
hypotheses  (5.1)  and  (5.2) 
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Thus  all  quantities  on  the  right  hand  side  of  the  inequality 

(5.15)  are  bounded,  and  hence  II  Ac  *  -  c  *||  is  bounded.  Therefore  a 

11  o  o 

finite  ’a*  exists,  for  A  :  0  <  A  <  1  ,  satisfying  the  fixed  point 

condition  given  by.  (5.12). 

» 

Now  consider  the  contraction  condition  which  states: 


|| Ac1*  -  Ac2*||  <_  A||Cl*  -  c2*||  (5.16) 

for  every  cj*>  c2*  £  S  ,  and  some  A  :  0  <  A  <  1  . 

Choose  some  cj_*>  c2*  e  S  .  Let  e^*,  e2*  t*ie  samPled 
error  functions  corresponding  to  c^*,  c2*  •  For  a  given  r*  it 
follows  from  (5.8)  that: 


p  *  _  r  sfe  s  p  *  —  a  * 

C1  C2  11  11  1  2 


(5.17) 


Also  from  (5.10): 


Ac^*  -  Ac2* 


B* 


(f>(e.*) 

- — - e  * 

K  1 


-  B* 


■<Ke2*) 


K 


-  e  * 
2 


<  B* 


<P  (ex*)  -  <j)(e2*) 


K 


-  (ex*  -  e2*) 


=  H* 


( e -  e2*) 


cf)(e1*)  -  ( He 2*) 


K  (ex*  -  e2*) 


-  1 


<_  H*||  e1*  -  e2* 


^(e^)  -  <j>(e2*) 


K  (ex*  -  e2*) 


-  1 


(5.18) 


From  (5.17)  and  (5.18)  it  follows  that  the  contraction  condition  will 


be  satisfied  if: 


. 


' 
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H*||e1* 


<l>(e1*)  -  c t>(e2*) 
K  ( e x*  -  e2*) 


1 


i .  e . 


4>  (e1*)  -  <j>(e2*) 

K  (ex*  -  e2*) 


(5.19) 


(5.20) 


i.e. 


i.e. 


A  :  *(el*}  -  *<e2*>  <2l. 

H*  -  K  (e  *  -  e2*)  1  -  H* 

AK  .  *<*!*>  -  ^e2*}  ,  AJK 
H*  -  e1*  -  e2*  -  +  H* 


(5.21) 


(5.22) 


Now  by  the  mean  value  theorem: 


<Ke  *)  -  <Ke2*) 

p  *  _  p  * 

1  2 


d<P(e) 

de 


e=e  *,  e  *  e  S 
m  m 


(5.23) 


Hence  (5.22)  will  be  satisfied  if: 


_  AK  d<Ke)  <  K 
H*  -  de  -  H* 


(5.24) 


Since  by  hypothesis  <  K2  ,  and  A  :  0  <  A  <  1  ,  it  follows 

from  (5.24)  that  the  contraction  condition  will  be  satisfied  if: 


K  -  f*  <  K1  and  K2  <  K  +  ~  .  (5.25) 

Hence  condition  (5.25)  alone  is  sufficient  for  the  mapping  A  to  have 
a  fixed  point  in  F  for  every  bounded  input  r  .  In  other  words 
system  N  is  BIBO  stable  if  (5.25)  is  satisfied.  This  completes 
the  proof  of  the  theorem. 
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5 . 4  Numerical  Example  1. 

To  illustrate  the  application  of  the  criterion,  consider  the 
nonlinear  integral-control  sampled-data  system  show  in  Fig.  5.2. 


c* 

c 


Fig.  5.2  Nonlinear  Integral-Control  Sampled-Data  System 


Let  4> (e )  be  replaced  by  some  linear  gain  K  ,  for  which  the 
resulting  linear  system  L  is  stable.  The  problem  is  to  pick  the 
optimum  value  of  K  which  gives  the  widest  possible  bounds  for  the 
slope  of  the  nonlinearity  for  the  BIBO  stability  of  the  system.  To 
do  this  the  behavior  of  H*  when  K  is  varied  is  to  be  studied.  The 
procedure  is  as  follows: 


The  Laplace  transform  of  the  total  linear  part  is  given  by: 


G(s)  - 


1  -  e 


-Ts 


s  (s  +  1) 


(5.26) 


The  corresponding  Z-transform  is: 


G(z)  = 


1  -  e 


-T 


z  -  1  -T 

z  -  e 


(5.27) 


% 

. 
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from  which  it  can  be  determined: 


H(z 

n  _  KG(z) 

'  1  +  KG(z) 

K[(T+e"T-l)z  +  l-e“T-Te_T] 

2  -T  -T  -T  -T  -T 

z  +  [(T+e  -l)K-l-e  ]z  +  K(l-e  -Te  )+e 

(5.28) 

Let  zv 

z 2  be  roots  of  the  characteristic  equation: 

2  -T  -T  -T  -T  -T 

z  +  [(T+e  -1)K- 1-e  ]z  +  K(l-e  -Te  )+e  =0  . 

(5.29) 

Since  L 

is  stable  K  should  be  such  that  |z^|  <  1  ,  and 

j  z  £  |  <  1  . 

In  order 

to  find  the  range  of  values  of  K  ,  for  which  |z^|, 

1 z2 i  are 

less  than  1  ,  we  may  use  the  following  criterion  applicable  to  a 

quadratic  polynomial  p(z)  ,  with  real  coefficients  and  the  coefficient 

2  25 

of  z  unity  : 

The  necessary  and  sufficient  conditions  for  p(z)  to  have 
no  zeros  outside  the  unit  circle  are: 

| p (0)  |  <  1  p(l)  >  0  p(-l)  >  0  .  (5.30) 


The  three  conditions  of  (5.30)  applied  to  the  polynomial  on 
the  left  hand  side  of  equation  (5.29)  lead  to  the  relations: 

|K(l-e"T-Te"T)+e“T|  <  1  (5.31) 

T(l-e_T)K  >  0  (5.32) 

2 (l+e~T) 

T(l+e~T)-2(l-e_T) 


K  < 


(5.33) 
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From  (5.32)  it  is  obvious  that  K  has  to  be  positive.  From  (5.31)  it 
follows  that: 


or 


— T  — T  — T 

-1  <  K( 1-e  -Te  )+e  <  1 


(5.34) 


1+e 


-T 


-T  -T 
1-e  -Te 


<  K  < 


1-e 


-T 


-T  -T 
1-e  -Te 


(5.35) 


Since  (5.32)  has  to  be  satisfied  (5.31)  reduces  to  the  requirement: 

1-  -T 

°  <  K  <  — br—yF  •  (5.36) 

1-e  -Te 

Hence  for  stability  K  should  be  such  that  (5.33)  and  (5.36)  are  satisfied. 
By  giving  various  values  for  T  it  can  be  seen  that  for  small  values  of 
T  ,  (5.36)  gives  the  upper  bound  for  K  ,  whereas  for  larger  values  of 

T  ,  (5.33)  gives  the  upper  bound  for  K  .  We  will  confine  our  discussion 

to  the  case  when  T  is  small,  i.e.  of  the  order  of  1  or  2  seconds. 

In  this  case  the  bound  on  K  for  stability  will  be  given  by  inequality 
(5.36).  For  this  range  of  values  of  K  ,  the  roots  z^,  z^  of  the 
characteristic  equation  (5.29)  may  be  real  and  unequal,  real  and  equal 
or  complex  conjugate.  By  equating  the  discriminant  of  (5.29)  to  zero, 
and  solving  for  K  we  can  determine  the  value  of  K  for  which 
z^  =  z^  .  The  value  of  K  thus  obtained  turns  out  to  be: 

~T 

K  =  -  •  ■  -  .  (5.37) 

(l+T-e“T)+2  /t  ( l-e_T) 

It  can  be  easily  verified  that  for  values  of  K  in  the  range  0  to  that 
given  by  equation  (5.37)  z^,  z^  are  real,  unequal  and  positive,  whereas 
for  values  of  K  greater  than  that  given  by  equation  (5.37)  upto  the 


} 
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limiting  value  for  stability  given  by  (5,36)  z^,  z^  are  complex  conjugate, 


Thus  we  can  say  z^,  z^  are: 


(1)  Real,  unequal  and  positive  if: 


0  <  K  < 


1-e 


-T 


(1+T-e  T)+2  A  (1-e  T) 


(5.38) 


(2)  Real  and  equal  if: 


K  = 


1-e 


-T 


(1+T-e  T)+2  v4(l-e  T) 


(5.39) 


(3)  Complex  conjugate  if: 


1-e 


-T 


<  K  < 


1-e 


-T 


(1+T-e  )+2 


-1*- 


.  -T  m  -T 
1-e  -Te 


(5.40) 


Consider  case  (1).  Let  z^  >  z^  .  Inverse  transforming 
H(z)  using  the  residue  theorem: 


(  0  ,  for  k  =  0  . 


h(kT)  = 


K  r  /m.  -T  1W  k  k.  .  -T  -T.  ,  k-1  k-ls  ,  . 

- [ (T+e  -1) (z  -z0  )+(l-e  -Te  )(zn  -z0  )]  .  (5.41) 

zj~z2  x  z  1  Z 

for  k  >  0  . 


Hence 


H*  =  l  | h (kT) |  =  l  h (kT)  . 
k=o  k=l 


(5.42) 


Note  that  the  absolute  value  sign  can  be  dispensed  with  in  summing  the 
series,  as  z^,  z ^  are  real  and  positive,  and  z^  >  ,  which  ensures 
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that  each  term  of  the  series  is  positive  for  all  k  =  1,  2,  .  .  . ,  00  . 
Moreover,  each  of  the  geometric  series  can  be  separately  summed  and  are 
convergent  because  of  the  fact  z  ,  are  less  than  1  .  Performing 

the  summation  we  get: 


l  h(kT) 
k=l 


K 


Z1"Z2 


'(T+e  T-l)z. 
1-z, 


(T+e  -X)f2  +  _l-Te-T-e_T  l-Te-T-e_T 


1-z , 


1-z. 


1-Z, 


K 


Z1"Z2 


(T+e  T-l)(2l-z2)  (1-e  T-Te“T)  (z^) 


(1-z^  (l-z2) 


( I- z -^)  ( 1—  ^2 ) 


KT(l-e  T) 

(l-z1-z2+z1z2) 


(5.43) 


Since  z^,  z2  are  the  roots  of  the  equation  (5029)  we  have: 


(Z;L+z2)  =  K(T+e  T-l)-e  T-1 


(5.44) 


and 


— T  -T  -T 
z  z  =  K(l-e  -Te  )+e 


(5.45) 


Substituting  (5.44)  and  (5.45)  in  (5.43): 

*  =  _ KT(l-e~T) _ 

ii  “  m  m  m  rn 

l+K(T+e  -l)-e  -l+K(l-e  -Te  )+e 

KT(l-e~T) 

KT(l-e"T) 


=  1 


(5.46) 


Now  consider  case  (2).  Here  z^  =  z2  =  zq  (say).  Application 


of  the  reside  theorem  in  this  special  case  yields: 
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h  (kT)  =  < 


0  ,  for  k  =  0 


K[ (T+e  T-l)kz  k  1+( 1-e  T-Te  T)(k-l)z  k  2]  for  k  > 
o  o 


.  ‘ .  H*  =  l  | h(kT) 
k=o 


-  K(T+e  T-l)  l  kz  K  X+K(l-e  X-Te  X)  £  (k-l)z 
k=l  °  k=l  ° 


k-1 


-T  m  -T, 


k-2 


Note  that 


l  kz_k  1  and  l  (k-l)z  k  2 
k=l  k=l  ° 


are  the  first  derivatives  of  the  geometric  series 


00 

V  k  , 

/  z  and 


k-1 


V  k-l 
)  z 
L  o 


k-1 


respectively,  and  hence  their  sums  are 


4-  br-]\  and  4" 

dz  1-z  dz  1-z 

z=z  ’ z=z 

o  o 


respectively.  Performing  the  differentiation  we  get: 


l  kz  k  1 
k=l  0 


(1-zo) 


and 


I  (k-1): 


k-2 


k-1 


(i-z  ) 
o 


2  ‘ 


(5.47) 
0  . 


(5.48) 


(5.49) 


(5.50) 


Substituting  (5.49)  and  (5.50)  in  equation  (5.48): 
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H*  = 


K(T+e  T-l)+K(l-e  T-Te  T) 

(1"Zo)2 


KT(l-e  T) 

l-2z  +z 
o  o 


(5.51) 


Since  z  =  z..  =  z~  we  have: 
o  I  2 


■2z  =  -(z.+zj  =  K(T+e  T-l)-e  T-1 

o  12 


(5.52) 


and 


2  w,  “T  m  -T.  ,  -T 

z  =  z..z0  =  K(l-e  -Te  )+e 
o  12 


(5.53) 


Substituting  (5.52)  and  (5.53)  in  (5.51)  we  get  as  before: 


H*  =  1 


(5.54) 


As  for  case  (3),  h(kT)  can  be  found  as  in  case  (1);  but 

it  is  difficult  to  get  an  expression  for  H*  in  a  closed  form.  However, 

an  estimate  for  H*  can  be  found  using  equation  (2.38).  One  may  also 

conjecture  that  H*  increases  from  1  to  ^  as  the  complex  roots  move 

from  z  to  the  circumference  of  the  unit  circle, 
o 


The  above  analysis  shows  that  as  long  as  K  is  such  that 
z^y  z 2  are  real,  H*  =  1  ,  and  the  permissible  bounds  for 


are  given  by: 


0  <  iilil  <  2K 

de 


(5.55) 


The  best  bounds  are  got  when  z^  =  z^  =  zq  ;  i.e.  when 


K  = 


1-e 


-T 


(5.56) 


(1+T-e  T)+2  A (1-e  T) 
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These  bounds  for  various  values  of  T  are  shown  in  Table  1. 

TABLE  1 

BEST  BOUNDS  FOR  d<^6^  FOR  BIBO  STABILITY 

de 

USING  CONTRACTION  MAPPING 


Sampling  period  T  (in  seconds)  12345 

Bounds  for  d-|(e)  :  0  <  <  0.42  0.31  0.26  0.22  0.19 

de  de 


5. 5  Comparison  with  some  Existing  Results: 

In  this  section  we  will  compare  the  results  obtained  above  with 

those  obtained  by  the  application  of  a  frequency  domain  criterion, 

3 

originally  proposed  by  Jury  and  Lee  for  the  absolute  asymptotic 
stability  of  a  class  of  nonlinear  discrete-data  systems,  but  later  shown 
to  be  sufficient  for  the  absolute  BIBO  stability  as  well  by  Iwens  and 
Bergen'7.  The  type  of  system  considered  by  them  has  the  same  configuration 
as  that  represented  in  Fig.  5.1(a),  but  the  constraints  on  the  nonlinear 
gain  function  <p(e)  are: 

<K0)  =  0  (5.57) 

0  <  <  K  (5.58) 

dcj)  (e) 
de 


<  K’ 


(5.59) 


5 
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It  may  be  mentioned  that  (5.59)  implies  that: 

K  <_  K'  .  (5.60) 

The  criterion  states  that  a  sufficient  condition  for  the  absolute 
asymptotic  as  well  as  BIBO  stability  of  the  type  of  system  considered 
is  that  the  following  relationship  is  satisfied  on  the  unit  circle  for 
some  non-negative  q  : 

Re  G  (z)  [  1+q  ( z-1)  ]  +  j  ^  |  (z-l)G(z)  |2  >_  0  .  (5.61) 

In  the  above  G(z)  is  the  Z-transfer  function  of  the  linear  plant  in  the 
forward  path  of  the  system. 

3 

Jury  and  Lee  in  their  paper  consider  the  same  system  that 
has  been  discussed  in  section  5.4  above,  as  an  example  to  illustrate  the 
application  of  their  criterion.  By  drawing  the  loci  of  G(z)  , 

(z-l)G(z)  and  (z-l)G(z)  -  ——  |(z-l)G(z)  | 2  for  z  -  exp  (jCDT)  and 
several  values  of  K'  ,  they  deduce  that  the  stability  criterion  for 
this  particular  example  can  be  reduced  to  the  following  inequality: 

-(§  +  D  +  £  +  q(T  -  Y”  T2)  >_  0  .  (5.62) 

Then  they  assert  that  if  K  and  K'  are  restricted  to: 

K  <.  K'  <  |  (5.63) 

there  always  exists  a  non-negative  q  for  which  (5.62)  is  satisfied. 

Thus  the  condition  for  stability  (both  asymptotic  as  well  as  BIBO) 
of  the  system  reduces  to  the  requirement  that  the  absolute  value  of  the 


i 

1 
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slope  of  the  nonlinearity  is  less  than  —  .  Table  2  gives  the  bounds 
on  the  slope  of  the  nonlinearity  for  stability  for  various  values  of 
T  .  For  comparison,  the  corresponding  bounds  on  the  slope  obtained  by 
the  contraction  mapping  technique  are  also  given  in  the  same  table. 


TABLE  2 

BEST  BOUNDS  FOR  d-|(e-  FOR  BIBO  STABILITY 

de 

(JURY’S  METHOD  VS  CONTRACTION  MAPPING) 


Sampling  period  T 

(in 

seconds)  1 

2  3  4  5 

Stability  theorem 

of  Jury  and  Lee 

|  d<J>(e)  |  r  o 

1  0.67  0.50  0.40 

'  de  ' 

Stability  via 

0 

<  dt(e)  <  0.42 

0.31  0.26  0.22  0.19 

contraction  mapping 

de 

It  must  be  conceded  that  the  results  obtained  via  contraction 

mapping  are  quite  conservative  compared  with  those  of  Jury  and  Lee. 

3 

In  fact,  Jury  and  Lee  state  in  a  foot  note  to  their  paper  ,  that 
Professor  Y.Z.  Tsypkin  has  indicated  that  condition  (5.59)  may  be 
replaced  by  the  more  general  condition: 

K 


—  00  < 


(5.63) 


} 
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which  gives  still  better  bounds  for  K'  for  stability.  However,  as 
far  as  the  contraction  mapping  method  is  concerned,  the  important  point 
to  be  noted  is  that  within  the  bounds  given  by  the  method,  the  actual 
solution  of  the  system  can  be  generated  by  a  process  of  successive 
approximations,  starting  from  the  solution  of  the  linear  system  as  a 
first  approximation,  and  using  equation  (5.10)  recursively;  and  the 
convergence  of  the  process  is  assured.  Jury’s  method  does  not  possess 
this  advantage. 


CHAPTER  VI 


GENERATION  OF  SYSTEM  SOLUTION 


6 . 1  Introduction. 

In  this  chapter  we  will  discuss  the  generation  of  the  solution 
of  a  nonlinear  discrete-data  system,  satisfying  the  contraction  and  fixed 
point  conditions,  by  a  process  of  successive  approximations.  It  will  be 
shown  that  the  process  results  in  a  discrete  Volterra  type  series,  whose 
convergence  is  assured  in  the  region  in  which  the  contraction  and  fixed 
point  conditions  hold. 

6 . 2  Derivation  of  the  Solution. 

To  illustrate  the  derivation  of  the  solution  consider  the 
discrete-data  system  with  a  single  cubic  type  of  nonlinearity  in  the 
feedback  path.  The  general  configuration  of  the  system  is  shown  in 


Fig.  6.1. 
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Fig.  6.1  Discrete-Data  System  with  a  Cubic  Nonlinearity 

in  the  Feedback  Path 
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It  will  be  readily  recognized  that  the  system  discussed  in 
example  1,  chapter  4.2  belongs  to  this  type.  For  convenience  the 
sampling  time  will  be  normalized,  so  that  for  the  rest  of  the 
derivation  it  is  understood  that  T  =  1  ,  The  functional  equation  for 
c*  for  the  system  is: 

c*  =  B*r*  -  B*Kc*3  .  (6.1) 

For  numerical  evaluation,  it  is  more  convenient  to  rewrite  the  above 
equation  for  the  output  at  the  end  of  the  nt31  sampling  period  thus: 

c(n)  =  Br(n)  -  BK[c(n)]3  (6.2) 

where  B  represents  the  convolution  summation  operator  for  that 
sampling  instant.  B  is  defined  by  the  relation: 

n 

Bx(n)  =  £  h(n-k)x(k)  (6.3) 

k=o 

where  h(k)  represents  the  impulse  response  of  the  linearized  system 

th 

at  the  end  of  the  k  sampling  period.  If  equation  (6.2)  is  solved 
for  c(n)  we  can  determine  the  sampled  outputs  at  the  various  sampling 
instants,  by  substituting  n  =  1,  2,  etc.  in  the  solution. 

In  deriving  the  solution  we  will  be  making  the  assumption 
that  all  initial  conditions  are  zero,  and  the  linearized  system 
corresponding  to  the  given  system  is  stable.  This  means  that  for  the 
linear  case,  i.e.  when  K  =  0  ,  the  solution  is: 

n 

c (n)  =  l  h (n-k) r (k) 
k=o 


(6.4) 
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Since  r(n)  is  such  that  the  contraction  and  fixed  point 

conditions  are  satisfied,  we  can  generate  the  solution  of  the  system 

by  a  process  of  successive  approximations,  starting  from  the  solution 

of  the  linearized  system  as  a  first  approximation,  and  using  equation 

(6.2)  recursively.  Let  the  output  of  the  linearized  system  at  the  end 
til 

of  the  n  sampling  period  be  denoted  by  cQ(n)  •  Then: 

n 

c  (n)  =  l  h (n-k) r (k)  (6.5) 

o  u 

k=o 


will  be  our  first  approximation  for  the  solution  of  (6.2).  Substituting 
(6.5)  on  the  right  hand  side  of  (6.2),  we  obtain  the  second  approximation 
for  c(n)  ,  viz.  c^(n)  . 


c1(n) 


n 

l  h(n-k) r(k) 
k=o 


-  K 


n 

I 

k=o 


h (n-k) 


k 

V 

L 

i_k^=o 


13 


h(k-k1)  r(kx) 


(6.6) 


The  last  term  in  (6.6)  can  be  written  in  the  following  way: 


n 


r  k 


K  l  h (n-k) 
k=o 


l  h(k-k1)r(k1) 


k^=o 


n 


K  l  h (n-k) 
k=o 


k  k  k 

l  h(k-k1)r(k1)  [  h(k-k2)r(k2)  l  h(k-k3)r(k3) 

Lk^=o  k2=o  ^3=° 


k  k  k  k 

K  l  l  l  l  h(n-k)h(k-k  )h(k-k  )h(k-k  )r(k  )r(k  )r(k„) 

k=o  k  =o  k2=o  k3=o 


n  n  n 

=  K  l  l  l 

k  5=0  k2=o  k3=o 


n 


l  h(n-k)h(k-k^)h(k-k2)h(k-k3) 
k  =  max 

L  k!»k2»k3 


r(k1)r(k2)r(k3) 


(6.7) 
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The  last  expression  in  (6.7)  is  obtained  by  rearranging  the  order  of 
summation. 

Let  us  now  define: 

n 

h^  =  l,  h(n-k)h(k-k^)h(k-k2)h(k-k3)  .  (6.8) 

k  =  max 

kl*k2’k3 

Then  (6.7)  becomes: 


n 


K  I  h(n-k) 


L 

k=o 


l.  k^=o 


l  h(k-k1)r(k1) 


n  n  n 

=  K  l  l  l  h3r(k1)r(k2)r(k3) 


(6.9) 


k^=o  k2=o  k3=o 


and  substituting  (6.9)  into  (6.6)  we  obtain: 


n  n  n  n 

c, (n)  =  l  h(n-k) r(k)  -  K  l  l  l  h  r (k1 ) r(k?) r (k  )  .  (6.10) 

k=o  k^=o  k2=o  k3=o 


If  we  desire  higher  order  approximations  we  must  then 
substitute  equation  (6.10)  into  the  right  hand  side  of  (6.2).  We  will 
then  obtain  the  solution  of  equation  (6.2)  in  the  form  of  a  discrete 
Volterra  type  series  as  follows: 


n  n  n  n 

c (n)  =  l  h(n-k) r(k)  -  K  l  [  l  h  r(k  )r(k  )r(k  ) 

k=o  k1=o  k2=o  k3=o 

+ 


+ 


c  •  • 


0  0  0 


(6.11) 
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The  above  series  is  analogous  to  that  obtained  by  Barrett 
for  a  continuous  system  with  a  single  cubic  nonlinearity.  The  system 
discussed  by  him  is  described  by  the  following  differential  equation: 

L(^)y(t)  +  ey3(t)  =  x(t)  (6.12) 

where  e  is  a  constant  and 

L(~)  =  L  (p)  =  pm+alPm_1+...+am  (6.13) 

a^,  a^,  .  ..,  a^  being  constants.  x(t)  is  the  continuous  input  and 
y(t)  the  output.  Letting  h(t)  be  the  impulse  response  of  the 
linear  part  of  the  differential  equation  (6.12),  and  assuming  that 
for  t  <  0  ,  x(t)  =  y(t)  =  0  ,  (6.12)  can  be  converted  to  a  nonlinear 

integral  equation: 

t  o  t 

y (t)  +  /  £h(t-u)y  (u)du  =  J  h(t-u)x(u)du  .  (6.14) 

o  o 

If  we  denote  the  convolution  operation  by  the  symbol  3  ,  the  above 
equation  can  be  written  as: 

y ( t)  +  3ay 3 (t)  =  3x(t)  (6.15) 

or  alternately: 

y ( t)  =  3x(t)  -  3ey3(t)  .  (6.16) 

Comparing  equation  (6,16)  with  equation  (6.1)  or  (6.2)  it 
will  be  readily  seen  that  they  are  analogous  to  each  other.  Proceeding 
by  the  method  of  successive  approximations,  Barrett  gets  the  Volterra 
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series  solution  for  the  equation  (6.14)  as: 


t  t  t  t 

Y (t)  =  [  h(t-u)x(u)du  -  /  /  /  £h3x(u1)x(u2)x(u3)du1du2du3 

o  o  o  o 


+  .  ..  ,  (6.17) 


where 

t 

h 3  =  /  h(t-u)h(u-u^)h(u-u2)h(u-u3)du  .  (6.18) 

o 


NOTE:  In  the  solution  given  by  Barrett  the  limits  of  integration  have 

been  shown  as  from  -  00  to  00  ;  however,  since  x(u)  =  0  for  u  <  0  , 
and  h(t-u)  =  0  for  u  >  t  ,  it  is  enough  if  the  integration  is 
performed  over  the  finite  range  0  to  t  . 

Comparing  (6.17)  and  (6.11)  the  analogy  between  the  continuous 
Volterra  series  solution  and  the  discrete  Volterra  series  solution  will 
be  evident.  Moreover,  the  contraction  mapping  principle  shows  that 
the  series  represented  by  the  equation  (6.11)  exists,  is  unique  and  is 
convergent  in  the  region  in  which  the  contraction  and  fixed  point 
conditions  hold. 

The  validity  of  the  Volterra  series  representation  of  the 

solution  of  a  discrete-data  system  can  be  justified  on  the  same  basis 

27 

as  has  been  done  for  the  continuous-data  systems  by  Christensen 
The  following  extract  from  reference  27  will  be  relevant  to  our 
discussion: 

"When  can  it  be  inferred  that  a  Volterra  series  exists  and 
is  convergent  for  a  given  system  if  the  contraction  mapping  principle 
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shows  that  system  to  exhibit  bounded  input  bounded-output  stability? 
The  answer  to  this  question  is  that  the  given  system  in  addition  to 
having  a  unique  solution  must  also  be  analytic. 

To  be  more  specific,  consider  the  usual  form  of  the  Volterra 

series : 


00  00  00 

( t )  *  /  h1(u)y(t-u)du  +  /  /  h2  (u-^u^y  (t-u^ 

—  00  _oo  _00 

y (t-U2)du^du2  +  ...  + 

00  00 

/  ...  j  h^(u^,...,  u_^  )y  (t-u^) .  .  . 

—  00  —00 

y (t-u_^) du^ . .  . du^  (6.19) 


which  gives  an  output  x^(t)  for  a  given  input  y(t)  to  the  system  in 

28 

question.  Brilliant  showed  that  if: 

00  00 

||  h^  ||  =  /  ...  /  Ih^Cu^,  . . .  .u^)  |  du^.  .  .du^  oo  (6.20) 

—  00  —00 

for  all  i  such  that  (6.19)  can  be  cast  into  the  form: 

00 

XX1  l  IlhJlY1  (6-21) 

i=o 

where  the  right  side  of  (6.21)  is  a  convergent  series  then  the  given 
system  is  an  analytic  system." 

For  a  discrete-data  system,  the  discrete  Volterra  series 
analogous  to  (6.19)  would  be: 


o'  08 Lb  .  roiii'Io  yu  r  gniv'kr: 


’ 
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n  n  n 

c(n)  =*  l  h1(k)r(n-k)  +  l  l  h2(k1,k2)r(n-k;L)r(n-k2) 

k»o  k^*o  k2=o 


n  n 

+  .  . .  ■♦■  £  ...  I  h  (k, 

k^=o  k^o 


,ki)r(n-k1) . . ,r(n-k  )+. . . 

(6.22) 


Note  that  in  (6.22)  the  integrals  have  been  replaced  by  summations. 
Moreover,  since  we  are  making  the  assumption  r(n)  =  0  for  n  <  0  , 
and  r(n-k)  •  0  for  k  >  n  ,  the  summation  need  be  performed  only 
over  the  finite  range  0  to  n  . 

(6.22)  can  also  be  written  as: 


c(n) 


n  n  n 

l  h1(k)r(n-k)  +  £  £  h^n-k^n-k^r^)  r(k2)  + 

k»o  k^=o  k2=o 


n  n 

. ..  +  £  ...  £  h  (n-k..  , . 

k^»o  k^-o 


.n-k^) r(k^) . . . r(k^)  +  ... 

(6.23) 


Note  that  (6.23)  corresponds  to  the  series  (6.11)  generated  by  the 
process  of  iteration. 

By  analogy  with  the  continuous  case,  the  condition  for 
analyticity  of  the  discrete-data  system  may  be  stated  as  follows: 

If 

00  00 

|| h  II  -  I  ...  I  |h1(k1 . k  )|  <  -  (6.24) 

k^»o  k^o 

for  all  i  such  that  (6.22)  can  be  cast  into  the  form: 


■ 
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II  c (n)  ||  _<  l  ||h  ||  ||r(n)||  1  (6.25) 

i=o 

where  the  right  hand  side  of  (6.25)  is  a  convergent  series  then  the 
given  discrete-data  system  is  an  analytic  system. 

Consider  now  the  series  generated  in  (6.11).  It  is  readily 
recognized  that  this  series  can  be  cast  into  the  form  (6.25).  Thus  the 
system  considered  is  analytic.  As  for  the  uniqueness  of  the  solution, 
the  contraction  mapping  principle  shows  that  the  solution  generated  is 
unique  in  the  region  in  which  the  contraction  and  fixed  point  conditions 
hold.  Hence  it  can  be  stated  that  a  convergent  Volterra  series  solution 
exists  for  the  system  considered. 
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CHAPTER  VII 


CONCLUSION 


7.1.  Summary. 

In  this  thesis  we  have  treated  the  bounded- input  bounded-output 
stability  of  certain  classes  of  nonlinear  discrete-data  systems  using 
the  methods  of  functional  analysis,  in  particular,  the  contraction  mapping 
principle.  After  discussing  the  stability  of  linear  discrete-data  systems 
from  a  functional  point  of  view,  sufficient  conditions  for  the  BIBO 
stability  of  discrete-data  systems  with  a  single  polynomial  type  of 
nonlinearity  have  been  developed  using  the  contraction  mapping  technique. 
The  results  have  been  stated  in  the  form  of  two  theorems  applicable  to 
two  different  types  of  systems  with  nonlinearities  located  in  different 
parts  of  the  systems.  A  method  of  computing  the  bounds  on  the  system 
inputs  and  outputs  for  stability  in  a  special  case  of  the  nonlinearity 
has  been  presented,  together  with  numerical  examples  to  illustrate  the 
technique.  Another  type  of  system  investigated  is  the  one  with  a  slope 
restricted  nonlinearity.  A  sufficient  criterion  for  the  BIBO  stability 
of  such  a  system  has  been  formulated,  and  its  application  also  illustrated 
using  a  system  which  has  been  treated  by  Jury  and  Lee  by  a  frequency 
domain  method.  Finally  a  method  of  generating  the  solution  of  a  non¬ 
linear  discrete-data  system,  satisfying  the  contraction  and  fixed  point 
conditions,  by  a  process  of  successive  approximations  has  been  given. 

It  has  been  shown  that  the  process  results  in  a  discrete  Volterra  type 
series,  whose  convergence  is  assured  in  the  region  in  which  the 
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contraction  and  fixed  point  conditions  hold. 

7.2  Conments  and  Suggestions. 

Since  the  contraction  condition  is  a  strong  requirement,  the 
stability  bounds  obtainable  by  this  method  are  bound  to  be  on  the 
conservative  side.  The  frequency  domain  method  is  superior  in  this 
respect,  as  is  observable  in  the  case  of  the  example  discussed  under 
slope  restricted  nonlinearity.  But  the  most  important  advantage  of  the 
contraction  mapping  method  lies  in  the  fact  that  once  the  stability 
region  is  found,  it  enables  one  to  generate  the  system  solution  by  a 
process  of  successive  approximations,  starting  from  an  initial  estimate 
within  the  region,  and  the  convergence  of  the  process  is  assured.  The 
frequency  domain  method  does  not  possess  this  advantage.  Moreover,  the 
frequency  domain  criteria  stipulate  that  the  nonlinearity  must  be 
entirely  contained  in  a  certain  sector;  they  thus  usually  refer  only  to 
global  stability.  The  contraction  mapping  method  can,  however,  be 
applied  to  study  local  stability,  particularly  when  the  nonlinearity  is 
explicitly  given  in  a  functional  form,  such  as  a  polynomial. 

Although  the  contraction  mapping  technique  presented  in  this 
thesis  deals  exclusively  with  a  single-input  single-output  system,  it 
is  believed  that  the  method  can  be  generalised  and  made  applicable  to 
multiple- input  multiple-output  systems.  It  is  also  believed  that  by 
choosing  suitable  norms,  it  should  be  possible  to  investigate  asymptotic 
stability,  as  well  as  stability  in  between  sampling  period.  These  are 
some  of  the  avenues  open  for  research  in  this  area. 


In  conclusion  the  author  wishes  to  state  that  the  main 
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contribution  of  this  thesis  in  the  field  of  Control  Engineering  is  the 
exploitation  of  some  of  the  basic  principles  of  Functional  Analysis  for 
the  investigation  of  the  BIBO  stability  of  a  class  of  nonlinear 
discrete-data  systems .  The  scope  of  the  thesis  is  perforce  confined  to 
theoretical  exposition,  though  the  author  realises  that  the  value  of  the 
investigations  could  have  been  enhanced  by  supplementing  it  with  some 
computer  simulation  results ;  this  could  perhaps  form  part  of  a  separate 
project  for  future  investigations. 

Some  attempt  has  been  made  to  bring  out  the  practical  importance 
of  BIBO  stability,  by  choosing  standard  practical  models  such  as  the 
integral  control  servo  to  illustrate  the  application  of  theory.  Since, 
for  practical  purposes,  the  mere  assurance  of  BIBO  stability  is  not 
enough,  an  attempt  has  also  been  made  to  obtain  an  estimate  for  the 
input  and  output  bounds  in  the  examples  discussed  in  Chapter  IV.  In  the 
present  state  of  the  development  of  the  nonlinear  theory,  it  is  indeed 
difficult  to  obtain  exact  bounds;  it  is  hoped  future  work  in  this  area 
will  help  to  resolve  this  difficulty. 

Finally  the  author  wishes  to  remark  that  while  the  practising 
engineer  is  apt  to  have  little  patience  with  theoretical  work,  it  is 
still  a  rewarding  field  as  even  the  purest  of  pure  mathematics  has 
the  habit  of  eventually  finding  its  way  into  application. 


aJiui.s  .  1  ' 
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